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CHAPTER I 
INTRODUCTION 

Computer simulation and radiation pattern analysis of 
fuselage-mounted airborne antennas at high frequencies are the object of 
this research. The primary goal of this investigation is to develop an 
analytic solution for predicting radiation patterns of airborne antennas 
in a more accurate and efficient manner. The radiation pattern analysis 
is based on the uniform geometrical theory of diffraction (UTD) which is 
a uniform extension of Keller's geometrical theory of diffraction (GTD) 
[1]. The UTD solution [2] for the radiation from sources on a convex 
surface is employed to compute the radiated field from the antenna 
mounted on an aircraft fuselage. The aircraft scattering analysis is 
based on the UTD solutions for wedge [3,4] and corner diffraction 
[6,21,29], The theoretical formulation of UTD is discussed in some 
detail in Chapter II. 
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One of the principle problems in the design of a reliable airborne 
antenna system is the location of the antennas on the aircraft structure 
in order to achieve the desired radiation coverage. If modern systems 
are to function properly, the antenna pattern must meet certain 
specifications. In fact, the antenna system performance is very much 
dependent upon the resulting antenna radiation patterns. Until 
recently, airborne antennas were located bv simolv choosina amona those 
locations allowed by the aerodynamacist which met the required 
specifications and scale model measurements then were used to evaluate 
the performance of the antennas in terms of its desired pattern. This 
approach of airborne antenna design requires a lot of engineering time 
and is also very expensive; furthermore, there are many antennas mounted 
on a single aircraft. Therefore, the need for an efficient analytic 
solution to compute airborne antenna radiation patterns is quite 
apparent. If these antennas can be located on the aircraft at the 
design stage using an analytic solution, then one can expect better 
performance in that more optimum locations and necessary structural 
changes can be anticipated. In addition, a future relocation or 
addition of antennas on an aircraft within its useful lifetime can also 
be easily accomplished through this approach. Once an optimum region is 
determined, the antenna can simply be flight tested to ascertain its 
actual performance. Not only can these analytic solutions be used to 
determine the best location, but they can, also, determine a more 
optimum antenna design for a given application. 
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Since it is a study of general -type aircraft, the 
aircraft is modelled in its most basic form. Previously the aircraft 
fuselage was simulated by an elliptic cylinder and the rest of the 
aircraft appendages were modelled by finite flat plates [7,8,9], 

However, the elliptic cylinder model could not predict the pattern close 
to the nose or tail sectors, where the deviation from the physical 
situation becomes very prominent. In fact, the fuselage has a dominant 
effect on the resulting radiation pattern of a flush-mounted airborne 
antenna. In order to overcome this deficiency, a composite prolate 
spheroid model was introduced [10,11]. The significance of this 
solution is that the spheroid model provides the proper polarization and 
curvature effects as opposed to the cylinder which models only one 
curvature. Note that the surface geometry dictates the polarization of 
the radiated field [2], However, the prolate spheroid representation of 
the fuselage is not general enough to satisfactorily approximate the 
wide variety of aircraft. The inadequacy asociated with the prolate 
spheroid results from Its circular cross-section. It has been shown in 
[8] that an elliptic cross-section is necessary to successfully simulate 
the wide variety of aircraft fuselage shapes. Therefore, the composite 
ellipsoid fuselage model used here becomes necessary. The development 
of the composite ellipsoid fuselage model should provide as much 
generality as needed to simulate the wide variety of aircraft found in 
practice. 
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In applying the UTD to antenna radiation problems involving curved 
surfaces, a major task is to determine the geodesic paths on the curved 
surface. An efficient numerical' technique to find the creeping wave 
geodesic paths on an ellipsoid [13,14] and the computer simulation 
technique of the whole aircraft structures are discussed in some detail 
in Chapter III, 

Based on these techniques, one applies a ray optics technique to 
determine the fields incident on the various scatterers. The fields 
diffracted are found using the UTD solutions in terms of rays which are 
summed with the geometrical optics terms. The rays from a given 
scatterer tend to interact with the other structures causing various 
higher-order terms. In this way one can trace out the various possible 
combinations of rays that interact between scatterers. Thus, one need 
only be concerned with the important structural scattering components 
which have dominant effects on the radiation pattern calculation, and 
neglect all other higher-order terms. The algorithms for the actual 
calculation of various UTD terms are discussed in Chapter IV. 

Using these algorithms, the Ellipsoid-Model Aircraft Code [15] has 
been developed to compute and superimpose the various UTD terms for the 
near or far field calculation. The validity of this analytic solution 
is verified in terms of comparisons with measured data for the wide 
variety of aircraft. The radiation pattern analysis is presented in 
Chapter V. 
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Since UTD is essentially a high frequency solution, the source and 
various scattering centers should be separated by at least a wavelength. 
In terms of the ellipsoid structure its semi -mi nor axes should be at 
least a wavelength in extent. In terms of the scattering from a plate 
this means that each plate should have edges at least a wavelength long. 
In addition, each antenna element should be at least a wavelength from 
all edges. In some cases even this requirement can be relaxed. The 
upper frequency limit is dependent on how well the theoretical model 
simulates the important details of the actual structure. 
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CHAPTER II 


THEORETICAL BACKGROUND 


A. INTRODUCTION 

The principal analytical tool employed in this report is the 
uniform geometrical theory of diffraction (IJTD). The UTD is basically a 
uniform extension of Keller's geometrical theory of diffraction (GTD) 

[1] but is valid within transition regions adjacent to shadow and 
reflection boundaries. The radiation from sources on convex surfaces 
and scattering from flat plates are analyzed using UTD in this 
chapter. The surface parameters of an ellipsoid which is used to 
simulate the aircraft fuselage are also introduced here. The curved 
surface and plates are assumed to be perfectly conducting, and the 
surrounding medium of the structure is free space. An exp (jut) time 
dependence is assumed and suppressed throughout this chapter. One 
is referred to references [16,34] for a more extensive treatment on the 
theory and application of UTD. 
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B. RADIATION FROM A CONVEX SURFACE 


The development of a uniform GTD (UTD) solution is presented in 
this section for the electromagnetic field radiated by an aperture in or 
a short monopole on a smooth perfectly conducting convex surface. For 
an aperture in a convex surface it is convenient to define an 
infinitesimal magnetic current moment dP rn (Q‘) at any point Q' in the 
aperture as 

dP m (Q') = E(Q')xn'da' (2.1) 


where 

E(Q') is the electric field, 

a 

n* is the outward unit normal vector to the surface, and 
da' is an incremental area at Q'. 

The tangential electric field in the aperture is assumed to be known. 

The dP m (Q') radiating in the presence of the perfectly conducting convex 
surface, which now covers the aperture as well, constitutes the 
equivalent source of the electric field dE m (P/Q') produced at any point 
P exterior to the surface. The total radiated electric field E m (P) is 
then found by integrating the Incremental field dE m ( P/Q ' ) over the total 
area A of the aperture. Thus, the radiated electric field is given by 

f n.( p ) - // <<P/Q'> • (2.2) 

A 
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Following the above development for the equivalent sources in the 
aperture radiation problem one may similarly define an infinitesimal 
electric current moment dP e (t') In dealing with the radiation by a 
monopole on a convex surface as 

dP e (Q') - KOdiV , (2 .3) 


where 


I(t') denotes the electric current distribution on the monopole 
and is assumed to be be known, and 

*' is the distance along its length measured from the base at Q'. 

It is also assumed that the monopole is a short thin wire whose total 
length h is much smaller than the principal surface radii of curvatures 
at Q'. The current moment dP e (t‘) radiating in the presence of the 
perfectly conducting surface constitutes the equivalent source of the 
electric field dE e (P/*') produced at P exterior to the surface. The 
total radiated electric field E e (P) then can be approximately 
calculated from a knowledge of only the field dE e (P/Q'), which is 
produced by the dP g (Q ' ) in the presence of the perfectly conducting 
convex surface, by simply replacing dP e (Q') by J dP e (z') cosCkt'cose 1 ) 
if P is in the lit region, or by dP a (i') if P is in the shadow 

region. It is noted that the cose^ inside the integral is defined by 

• A A A 

cose * n' • s, where s is the radiation direction in the lit zone from 
any point on the monopole, with 0 < i* < h. 
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According to geometrical optics, the space surrounding the source 

is divided into an illuminated and shadow regions by a plane tangent to 

the surface at Q'. This plane is referred to as a shadow boundary. The 

present formulation or ansatz of the uniform GTD solution leads to 

separate representations for the radiated field dE m (P/Q‘) in the shadow 

e 

and lit regions, respectively. However, these different representations 
will be shown to match exactly in polarization, amplitude, and phase at 
the shadow boundary. 

1. Shadow Region 

According to the generalized Fermat's principle, a ray emanating 
from a source, which is located on the surface, follows a geodesic path 
on the surface and continually sheds energy into the shadow region. 

Such a creeping wave mechanism is illustrated in Figure 2.1(a), from 
which it can be seen that a surface ray traverses from the source point 
at Q', follows a geodesic path TH) to the diffraction point Q and then 
propagates along the geodesic tangent at Q toward the observation point 
P $ which is located in the shadow region of the curved surface. As 
shown in Figures 2.1(b) and (c), n' is the unit surface normal vector, 
t' is the unit tangent vector of the geodesic path and the binormal 

AAA 

unit vector is given by b* * t* x n' at the source point Q'. Similar 
unit vectors are introduced along the geodesic path. Thus at the 

a A 

diffraction point 0, n is the unit surface normal vector, t is the unit 
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tangent vector of the geodesic and the bi normal unit vector is given by 

AAA 

b=txn. The following solution Is valid for torsional geodesic paths, 

A A 

where b'*b. 

The incremental field dE m at P s can be expressed in terms of the 

e 

field at a reference point P Q by [2] 


«m(%) ' 0E m (P ) 
e e 


d d 
p l p 2 


-jks # 


'VV (o 2 t5 o' 


+ 0[m 


-2 


m" 3 ] 


(2.4) 


where and are the principal radii of curvatures of the wavefront 
-2 -3 

and 0[m , m 1 are the higher order terms. 

It is seen that if the reference point P q is moved t0 the curved 

surface diffraction point Q, then Pj -► o, Pg ♦ p c , and s Q ♦ s. Note 

that p c is the surface diffracted ray caustic distance which can be 

found via differential geometry [17]. Since dE m (P s ) is independent of 

e 

the reference point P Q , such that 


"* Jl~ W * *>•<» •• 

* 0 ^ e e 

pJ+0 


(2.5) 


then. 
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dE (P ) ~ L (Q',0) —7 r e jkS 

m v s' t m vw , s ( p.+s ) e 

e e y c ' 


Furthermore, L m (Q\Q) can be related to the source strength dP 
e 

by 


L mtQ',0) ■ dP n (Q').yq',Q) 
e e 


"nr 

e 


where T m (Q',Q) is given by [2] 


-jk A . 

% 

A A 

a 4u L b * nTj (Q* )H 

+ t'bT 2 (q 

e -jkt 

1 % 

Pg(Q) 

\ 

1 dn(Q) 

Pg(Q') 

-jkZ 0 



■ 4tt l n ' nT 5(Q' 

1 )H + n'bT 

| 

! d* 

Pn(Q) 

e-jkt 


! 9 

J 

dn(Q) 

Pg(Q') 


1/6 


1/6 


( 2 . 6 ) 
at Q' 

(2.7) 

4<Q* )Hl 

( 2 . 8 ) 

(2.9) 


12 


The quantities T^Q') , T 6 (Q* ) are the torsion factors at Q' 

and are given in Table 2.1. It is noted that the effect of torsion is 
confined to the source location and depends on the direction of the 
surface ray with respect to the principal directions of the surface at 
the source point. In other words, it does not depend on the surface ray 
path away from the source. The uniform Fock functions are expressed as 


H = g(5) 


and 


-j 

S = m(Q') 


gU) 


with 


1 ” 

9(5) - /7 / 

«exp(-j2ir/3) 


dx 


exp(-jx£) 

w‘(t) 


( 2 . 10 ) 


( 2 . 11 ) 


( 2 . 12 ) 


and 


g(o - n l 


dx 


«exp(-j2ir/3) 


exp(-jxS) 
wJ(t) 


(2.13) 


which are known as the acoustic hard and soft Fock functions, 
respectively. The Fock type Airy function is given by 
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(2) The quantities E and G denote two of the three coefficients 
E, f, G that appear In the “first fundamental form" of 
Differential Geometry [17 ], 


























X 00 

w ? (t) = /; f dt.exp (Tt-t 3 /3) 

“exp(j2ir/3) 


(2.14) 


and # 2 ( 1 ) is the derivative of Wgfx) with respect to t. The Fock 
parameter £ for the shadow region is given by [2] 


Q m(t’) 

£ = / dt 1 n • \ 
Q' p g^ ' 


(2.15) 


with 


1/3 


m( t ' ) 


kPg(t') 


(2.16) 


Here Pg(x') is the surface radius of curvature along the ray path at t'. 
The wavenumber (k) is taken to be that of free space. The width of the 
surface ray tube at Q,dn(Q), is given by 

dn(Q) * P c dY . (2.17) 

Note the parameters Z 0 and t are defined as the free space wave impedance 
and geodesic arc length from Q’ to Q, respectively. 
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A A 

Combining Equations (2.6) - (2.17), the n and b directed components 

of dE m (P $ ) are given by 
e 


(a) dP m (Q* ) case: 


dE^Pj (dP m .b‘)He" jkt 


f H in 


Pg(Q’) 

l d * 0 

, 

i 

- fn\ 
MgVV/ 

1 A, i. 

1 

V 

r f ^ ic \ 
* * 7 


e" jks + 0[nf 2 ] 


(2.18) 


and 


<( p s> 




-jkt 


Pg(Q') 

W 


- 1/6 



1 

s(p+s) 
H L 


-jks _ 2 .3 

e + 0 [m , m ] 


(2.19) 


(b) dP e (Q') case: 


-jkZ 0 . -jkt 

Pg(Q') 

K 

> 

* 4 it dP e (Q')He J 

_P g (Q) _ 

/ d * 

J s(p c +s) 


e _jks + 0[m‘ 2 ] 


( 2 . 20 ) 
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and 


dE e< p s> 


-JkZ 0 

4tt 


dPe(Q')T 0 S e:jkt 


Pg(Q') 

p g (Q) 



J s ( p c +s ) 


e" jks + OCm- 2 ] 


( 2 . 21 ) 


where T 0 = T(Q')pg(Q') with T(Q') being the surface torsion at the 
source location (refer to Table 2.1). The ray divergence factor 
Aiity)/di|> is used to describe the amount of separation between adjacent 
geodesic paths as surface rays propagate around the surface as shown in 
Figure 2.1(b), in which d\j» 0 anc * di|» are the angles subtended by the 
surface ray tube at the launching Q' and diffraction Q points, 
respecti vely. 

2. Lit Region 

In the lit region, the radiated energy propagates directly from the 
source Q' to the field point P|_ along the incident ray of geometrical 
optics as shown in Figure 2.2(a). For convenience, a ray-fixed 
coordinate system is used to describe the polarization properties. The 

A A 

unit vector s is in the direction from Q' to P, and n' is again the unit 
surface normal vector at Q'. The incident plane is defined as the plane 

A A A 

containing the unit vectors s and n'. The unit tangent vector t' is the 

A 

projection of the unit vector s In the tangent plane such that 
s * n'cos© + t^sinQ . The tangent plane is defined as the plane 
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(b) 


Figure 2.2. Ray tube and ray coordinates for the lit region. 
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tangent to the curved surface at Q'. Hence, t is tangent to the 
curved surface at Q'. The binormal unit vector at Q' is defined by 

A I A I A A' 

b^ = t^ x n'. The unit vector n is in the incident plane, perpendicular 

A 

to s and directed outward from the curved surface. The bi normal unit 

AAA 

vector at is given by b = s x n. 

In the ray optical approximation, the field dE m ( Pj_) at the field 

e 

point P[_ may be expressed by 


dE (P, ) ~ dEfPj 
m v L' nr o 





( 2 . 22 ) 


Since Q' is the only caustic of the incident rays, the principal radii 

i i ~ 

of curvature and p 2 associated with the incident wavefront at P Q are 

identical, i.e., p| 8 pj s p 1 . Furthermore, dE m ( Pj_ ) is independent of 

•m fSr 0 

the reference point P Q . If P Q is chosen to be at Q' , it follows that 

11m ^ 
e e 

p 0 

p 1 ♦ o 

r o * s (2.23) 

should exist. Thus, L* can be related to dP m (Q ’ ) by [2] 

e e 
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(2.24) 


L* = dP ( Q ' ) • T* 
m m ^ ' m 

e e e 


Equations (2.22) - (2.24) are, then, combined to yield 


di m< p L> -<">') 
e e 


T* 

m 


»-jks 


°[m: 2 . 


* • • • ] 


where 


m 


jk A I A A I A A I A Al A 

f b * nA + v 5 + b A bc + t Jt nD ] 


and 


-jkZ 0 

“77" [ n ‘ nM + n * bN 1 


where A, B, C, D, M, and N are defined in Table 2.2. Note that 

A A 

is decoupled into n and b components as follows: 

(a) d P m (Q') case: 


<( P L) 


- j k 

* 7T [MP.VltH'+TjFcos# 1 ) ♦ (dP.-t* IT/cose 1 ] 
-jks 


e 


s 


+ 



and 


(2.25) 

(2.26) 

(2.27) 

dE >L> 


(2.28) 
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and 


k -jkZ 0 . e -jks 

dE e< p |J s — dP e^') sin9 T o F — 


+ 0[m“ 2 ] . 


(2.31) 


The functions H* and S* are defined by 

. -j d/3 

H* = g(5 4 ) e 


and 


S * = m(Q') e 


-J5?/3 


(2.32) 


(2.33) 


in which the hard and soft Fock function g and g have been defined 
previously in (2.12) and (2.13), respectively. The Fock parameter c 
for the lit region is given by 


= - m )l (Q , )cose 1 


(2.34) 


where 


m(Q') 

m t (Q,) = /. T 2 2 i, 1/3 
(l+T 0 cos 0 ) 


(2.35) 


i * a .• 

Here the angle 0 is defined by n'»s * cos0 as shown in Figure 2.2. 
Also, 

S* - H^cos© 1 


I+TqCOS 2 © 1 


(2.36) 

as defined in Table 2.2. The other parameters are the same as before. 
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3. Pattern Factor 


The solutions for short magnetic or electric dipoles have been 

given In Sections 1 and 2. One approach to analyze an extended aperture 

or linear antenna problem is to integrate the above solutions over the 

source distribution as discussed previously. This is an application of 

the superposition theorem, and one approximates the source distribution 

by an array of short magnetic (or electric) dipoles on the conducting 

surface. This is an accurate solution, however, it is rather tedious. 

A more efficient approach has been introduced in reference [18]. It is 

assumed that the dominant waveguide mode exists across the aperture 

(i.e., a cosine distribution along the B dimension and uniform along the 

A dimension). The slot can be axially oriented or circumferentially 

oriented as shown in Figure 2.3. Notice that any slot with arbitrary 

orientation can be resolved into these two orthogonal components. 

For the monopole case, a sinusoidal distribution with unit amplitude is 

assumed for the current density over the monopole of length L. This 

approach is to modify dP m (Q') such that 

e 


(a) in the shadow region: 


m 


* 2B 

^m ir 


cos dr (V* 1 )) 

/ kB * * 2 


sin (T P m *b ) 

(oTT a 

— P m .b‘ 

2 m 


(2.37) 
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Figure 2.3. A conducting cylinder upon which (a) an axial slot or 
(b) a circumferential slot is mounted. 


and 


P e = n' [l-cos(kL) ] 


(2.38) 


(b) in the lit region: 


. 2B 

D _ n 

.kB i * * 

cos (y sine 1 (P m *t‘)) 


,kA i * 

sin (y sine^P^b')) 

r - r_ 

tn m it 

,kB i « « 2 

l-(— sine 1 (P m *t')) 


kA * ~ 

T Sin9 


and 


P e = n' 


cos(kLn' •s)-cos(kL) 

^ * 2 
l-(n'»s) 


(2.39) 


(2.40) 
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^ ere * P m * un1 *t vector in the direction of magnetic current moment, 

A,B are dimensions of the slot, and L is the length of the monopole. It 
is noted that L is not to exceed a quarter wavelength for the solutions 
to be valid. 

C. SCATTERING FROM FLAT PLATES 

1. Diffraction by a Wedge 

The Uniform Geometrical Theory of Diffraction (UTD) developed by 
Kouyoumjian and Pathak [3,4] is sufficiently general to handle the 
three dimensional wedge structures. The UTD wedge diffraction 
formulations are based on the fundamental Geometrical Theory of 
Diffraction (GTD) which was originally developed by Keller and his 
associates at the Courant Institute of Mathematical Sciences [1]. The 
GTD is a ray optical technique and, therefore, it allows one to gain 
substantial physical insight into the significant physical mechanisms 
involved in the scattering of flat plates. Accordingly, one is able to 
determine the dominant diffraction mechanisms for a given geometry. 

The field radiated from a source at 0 and observed at P in the 
presence of a perfectly conducting wedge is shown in Figure 2.4. It is 
assumed that the source and observation points are sufficiently removed 
from the wedge surface so that the contributions from the surface ray 
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field can be neglected. Applying the generalized Fermat's principle, 
the distance along the ray path OQ E P is a minimum which in terms of edge 

A A A A 

diffraction this leads to s'»e = s»e. 

Kouyoumjian and Pathak have shown that the diffracted fields can be 
written in a compact way if these fields are written in terms of an 
edge-fixed coordinate system [19,20]. This coordinate system is 
centered at the diffraction point Q E , and this point is unique for a 
given source and observation points. The incident ray diffracts as a 

I 

cone of rays such that & Q = B Q . The orthogonal unit vectors associated 

A I A A 

with the edge-fixed coordinate system are defined as S ■ <(>• x s' and 

AAA A 

B = $ x s, where s' is the unit vector in the direction of incidence, 
o 

A 

and s is the unit vector in the direction of diffraction. 

The diffracted field Eq can be asymptotically expressed by [3] 

E d (s) - E 1 (Q e ) • D(s,s')A(s)e‘ jks (2.41) 


in which D(s,s') is the dyadic edge-diffraction coefficient as given by 

5 a -*C B o D s " *’* D h . (2.42) 

The diffraction fields can also be written in a matrix form as 


El<s) 

S 

■°s 0 


E «< 0 e) 

*!<») 
I_ _ 


1 

0 

1 

O 

l =■■ 


e i (5 e> 


A(s)e 


-jks 


(2.43) 
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Expressions for the scalar diffraction coefficients which are valid at 
all points away from the edge (again excluding * 0 or mr) are [3] 


»/4 


D s,h (L, * , *'* B o ) “ 2n/HJTsine 


r ir+6 


cot (^-)F(kLa + (n) 




+ cot (Tn-)F(kLa-(B-)) * 




cot (- 1 2^-)' : (kLa + (0 + )) 


,tt-B 


+ cot (-^p)F(kLa"(3 + )) 


(2.44) 


where the minus (-) sign applies to D s and the plus (+) sign for D^. 
Note that D s is referred to as the soft scalar diffraction coefficient 
for the acoustically soft (Dirichlet) boundary condition at the surface 
of the wedge, i.e.. 


( E 


) * 0 


Wedge 


(2.45) 


Dj, is referred to as the hard scalar diffraction coefficient for the 
acoustically hard (Neumann) boundary condition at the surface of the 
wedge, i.e.. 



) * 0 

Wedge 


(2.46) 
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where the 3/3n is the derivative along the normal direction to the 
boundary surface. The angle parameter 0 is given by 


0 + » $ * 


(2.47) 


and 


±. ±. 
a (6 ) 


2cos 


. 2 ^- 0 * , 
( 2 ) 


(2.48) 


where N 4 are the integers which most nearly satisfy the following 
equations: 

2irnN + - 0* = ir (2.49) 


and 


2irnN“ - 0~ = -u . (2.50) 

Note that the 0" terms are associated with the incident field and the 
0 + terms are associated with the reflected field as shown in Table 2.3. 
The wedge angle number (n) is given by 

WA 

n ■ 2 - — » (2.51) 

where WA is the wedge angle in radians. The wedge transition function 
which is basically a Fresnel integral is given by 
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TABLE 2.3 

COTANGENT FUNCTIONS FOR THE DYADIC DIFFRACTION COEFFICIENT 




The cotangent is singular when 

value of N at the boundary 

cot 



AT -0 

\ 2 rt } 

surface ^ » 0 is shadowed 

cot 


4> * 4>‘ + n, a SB 

AT -0 

l 2n I 

surface ^ * nn is shadowed 

cot 

( n + + <£') \ 

l 2 n ) 

6 = (2n — a RB 

reflection from surface <p » nn 

N* - 1 



4> * n - 4>\ a RB 

N~ -0 

tui 

\ 2 ft I 

reflection from surface ^ * 0 


F(x) = 2j/x eJ x e"J x2 dt . (2.52) 

The magnitude and phase of this transition function F(x) are shown in 
Figure 2.5. 

When x is small, F(x) is given by 

F(x) ~ [/irx - 2xeJ' ir / 4 - 2/3 x 2 e"^ 4 ] eJ^/ 4 + x ) (2.53) 

and when x is large, F(x) is given by 
j 3 1 15 1 75 1 

F(x) ~ ! + 27 - -j -Q^r + yg^- . (2.54) 
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Figure 2.5. Transition function. 
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PHASE ( DEGREES ) 


When 0.3 < x < 5.5, i.e., in between those regions, a linear 
interpolation scheme is used such that 

F(x) * F(x i ) + (F(x i+1 ) - F(x i )) (x 1+1 - Xl ) (2.55) 


F( x i+l)" F ( x i ) 

where the F(x,.), — — and x, can be tabulated. These 

v i+l" A i' ' 

solutions can be used to conveniently compute the transition function. 
The L is a distance parameter, which is dependent on the type of 
illumination, and is given by 


. 2 

s sin 8 

0 

for plane wave incidence 

ss' 


s+s' 

for cylindrical wave incidence, and 

2 

ss' sin 0 O 

for conical and spherical wave incidence 

s+s* 


The spreading factor A(s), which accounts for the spreading of the power 
along the diffracted ray, is given by 

plane, cylindrical and concical wave 
incidence, and 

spherical wave incidence 

(2.57) 
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2. Diffraction by a Corner 


A corner is formed by the intersection of a pair of finite straight 
edges. Figure 2.6 shows the corner diffraction geometry. The 
diffraction of energy from these edges is based on diffraction from an 
infinite straight edge. To compensate for the finiteness of the edges, 
a diffraction coefficient associated with the corners of the plate 
is needed. An empirical corner diffraction solution proposed by 
Burnside and Pathak [21] is based on the asymptotic evaluation of the 
radiation integral which employs the equivalent edge currents [6,31] 
that would exist in the absence of the corners. The corner diffraction 
term is then found by appropriately (but at present empirically) 
modifying the asymptotic result for the radiation integral which is 
characterized by a saddle point near an endpoint [28]. Even though 
this corner diffraction coefficient is still in its development stages, 
it has been shown to be very successful in predicting the fields 
diffracted by a corner for a number of plate structures. For this 
reason, it is discussed here and has been used to obtain some 
interesting results. 

The corner diffracted fields associated with one corner and one 
edge in the near field with spherical wave incidence are given by [21] 


X 


1 Z o 

s 


z c . 


M Y 0 





/ sin & c sin B 0C ^ ^ e‘ jks 

(cos8 oc -cose c ) F f kL c a ^ +6 oc" 8 c^ 4ns 


(2.58) 
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RECEIVER 

POINT 


Figure 2.6. Geometry for corner diffraction problem. 


where 


I 




c s (Q E )Y 0 

M 


E V«c>' 


W Z o 


/ 


■ST 

IT 


-jir/4 


_ e -j"/4 

C s (q E> * 27?SsinF 
h o 


F(kLa(8")) 

cosCg”) 


F { 


La(8~)/X 

kL c"(’ +B oc- B c> 


) 


T 


F(kLa(9*)) 

cos(j~) 


La(8 + )/X 

F kL a(ir+e -p )} 
c oc c' 


(2.60) 


The angle parameters 8 C and 8 0c are shown in Figure 2.6. The function 
F(x) was defined in the previous section, and 

a( 8) = 2cos 2 (8/2) 

6* » <j> ± 4.' 
s's" 

L a Is' + s") 

for spherical wave incidence. The function C s> h(Q) is a modified 
version of the ordinary edge diffraction coefficient for the half-plane 
case (n=2). The modification factor, which is given by 


(2.61) 


v 

s1n 2 S 0 and L c . — 



is a heuristic function that ensures that the diffraction coefficient 
will not change sign abruptly when it passes through the shadow 
boundaries of the edge, i.e., the corner diffracted field will ensure 
the continuity of the fields as the edge diffracted field shadow 
boundary is crossed. There is also a corner diffraction term associated 
with the other edge forming the corner and is found in a similar manner. 

D. ELLIPSOID SURFACE PARAMETERS 

In applying UTD to electromagnetic diffraction problems, some 

elementary knowledge of differential geometry of surfaces is necessary. 

The fuselage of an aircraft is modeled as a composite ellipsoid in order 

to better approximate general fuselage shapes. Using the ellipsoid 

* 

geometry shown in Figure 2.7, the surface is defined by 

*► A A 

R(e,<f>) = R(9 f <j))sin9cos4>x + R(0,<f>)sin9sin<j>y + R( 9 >( f>)co$0z 

(2.63) 

or, 

R(V e ,V r ) * a cosV e cosV p x + b cosVgSinV^ + c sinV g z (2.64) 


with 


cote 


tan V g = 


cos 2 <j> 


5 

si n $ 

~b* 


and tan V„ * 


asin<f> 


r bcos <|> 


(2.65) 
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Figure 2.7 


Geometry of an ellipsoid. 



The V p and V g parameters are introduced because of their convenience in 
analyzing elliptic geometries. 

Considering a ray which propagates along a geodesic path on the 
ellipsoid surface as shown in Figure 2.8, the three unit vectors t, n, 

A 

and b are, as defined earlier, the geodesic tangent, outward surface 
normal, and binormal at any point along'the geodesic path. The outward 

A 

surface unit normal (n) is obtained from 


n 



( 2 . 66 ) 


3R 

A A A 

= -a sinV e cosV p x - b sinV e sinV p y + c cosV z 
e 

and 

-*• 3R A A 

R V = Tv7 = " a cosV e sinV r x + b cosV e cosV p y 


where 



(2.67) 


( 2 . 68 ) 


Then, 


n * 


2 * 2 * * 
be cos V e cosV r x + ac cos V e sinV r y + ab cosV e sinV e z 

[b 2 c 2 cos 4 V e cos 2 V p + a z c 2 cos 4 V g sin 2 V p + a 2 b 2 cos 2 V g sin 2 V e ] 1/z 


cosV e cosV r „ cosV e sinV r A sinV e * 
" a x 4- b — ■ y + — c — z 

A 


(2.69) 
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where 


1/2 



cosV e cosV- 2 cosVp s1nV r 2 sinV- 2 

(— *i — L ) * ( — 1 — £ ) M-r 1 ) 


(2.70) 


The normal curvatures K n on the surface are evaluated by introducing the 
first and second fundamental forms of differential geometry [17] such 
that 

L dVp + 2M dV r dV e + Ndv| 

n E dV^ + 2F dV r dV e + GdV* ^ 2 * 71) 

^ ~ *v V • N t M * Rw w • N , N = Rw u • N 
r r r e ¥ e e 

(2.72) 

E - Ry • Rw , F * Rw • Rw » G - Rw • Rw 
r r r e e e 


A A 

and N * -n. It can be shown that 

R v v = -a cosV cosV x - b cosV sinV y - c sinV e z 
e e e 

n a A 

Ry y = a sinV e sinV r x - b sinV e cosV p y 
e r 

and 

■* * A 

R v v * -a cosV e cosV p x - b cosV e sinV r y 


(2.73) 
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After some algebraic manipulation, one obtains 

cos 2 V e 
L = “A 

M * 0 

1 

N " A 

E = a 2 cos 2 V g sin 2 V p + b 2 cos 2 V g cos 2 V p 

F = (a 2 - b 2 ) sin 2 V g sin 2 V p (2.74) 

G = a 2 sin 2 V g cos 2 V p + b 2 sin 2 V g sin 2 V p + c 2 cos 2 V g . 

A pair of orthogonal directions exists for which curvature, K, 
assumes maximum and minimum values, i.e., principal directions 

A 

A 

represented by two unit vectors Tj and t 2 . The two extreme values of K 
corresponding to the above directions are called principal curvatures 
denoted by and K 2 [26], The mean curvature is given by 


Kj + K 2 EN - 2MF + LG 
K M = 2 = 2(EG - F*) 

and the Gaussian curvature Is 

LN - M 2 
K g * K 2 * EG _ F Z 


(2.75) 
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Thus, the principal curvatures K* and K 2 are given by 


K 


1.2 


* K M * 



The two principal directions (tj, T 2 ) are given by 


(2.76) 


1 


1 K + a R V 
r e 


and 


(2.77) 


T 2 = To 


S R>. + 1 Rw 

r e 


where 


L - KjE 
a = K X F - M 


M - K 2 F 
* 8 = K 2 E - L 


(2.78) 


o 1/2 2 1/2 

* (E + 2aF + a G) , and y 2 « (BE + 26F + G) 


However, it is noticed that q ~ L/E and K 2 ~ N/G within the significant 
energy region for most practical cases. That indicates approximate 
values of the two principal curvatures, Kj ~ L/E and K 2 ~ N/G are good 
enough to be used for the geodesics on the ellipsoid surface for our 
radiation consideration. Thus, one obtains 
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K 


1 



cos 2 V 

e 

A[a 2 cos 2 V a sin 2 V„ + b 2 cos 2 V„ cos 2 V_] 
L e r e r J 

(A[a 2 sin 2 V p + b 2 cos 2 V r ]) 1 


and 


K 


2 


N 1 

® A[a 2 sin 2 V cos 2 V_ + b 2 sin 2 V sin 2 V + c 2 cos 2 V 1 
u e r e r e ' 


(2.79) 


(2.80) 


It is noticed that the principal radii of curvature are by 

* “ A A 

Rl = "jq and R 2 = "j^. If 6* denotes the angle between t and then 

A A A 

t = cose' + Tg sine 1 . From Euler's theorem, the normal curvature 
along the geodesic path is specified by 

Kg = KjC0S 2 e* + KgSin^' (2 

1 

with the radius of curvature, p , being . The torsion term (T ) 

5 ? 3 

introduced in Section (B) is given by 


T o = T * p g (2.82) 

where the surface torsion is given by 
sin 2 e' 

T a ~ ( k 2 " (2.83) 

with K^ and Kg being defined in Equations (2.79) and (2.80). 
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CHAPTER III 


FUSELAGE SIMULATION TECHNIQUE 


A. INTRODUCTION 

For airborne antennas mounted on the fuselage of an aircraft, the 
fuselage has been previously modeled as an elliptic cylinder [7,8,9] or 
a prolate spheroid [10,11]. However, the elliptic cylinder model could 
not predict the pattern close to the nose or tail sectors, where the 
deviation from the physical situation becomes very prominent, and the 
prolate spheroid model is not general enough to satisfactorily 
approximate the wide variety of aircraft cross sectional shapes. The 
inadequacy associated with the prolate spheroid results from its 
circular cross-section. It has been shown in reference [8] that an 
elliptic cross-section is necessary to successfully simulate the wide 
variety of aircraft fuselage shapes. In addition, the dominant 
structural effect in the elevation plane pattern is the profile of the 
aircraft fuselage. Note that the UTD solution in the lit region does 
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not depend on the surface parameters in that it is assumed the source is 
mounted on an infinite ground plane tangent to the surface at the source 
point. On the other hand, the transition and deep-shadow region 
solutions are modified due to their dependence upon the surface 
parameters. Therefore, the composite ellipsoid fuselage model becomes 
necessary and will be discussed in detail in Section III-B. In fact, 
the fuselage has a dominant effect on the resulting radiation pattern of 
a flush-mounted airborne antenna. 

In applying the UTD to antenna radiation problems involving curved 
surfaces, a major task is to determine the geodesic paths on the curved 
surface. Among the solutions for obtaining the goedesic paths of the 
three models, i.e., the cylindrical, prolate spheroidal and ellipsoidal 
models, the ellipsoid case is the most involved and complex one. This 
is to be expected because the equation describing an ellipsoid can 
degenerate into that of a cylinder of a prolate spheroid by using 
appropriate simplifications. 

According to the generalized Fermat's principle, a ray emanating 
from a source, which is located on the surface, follows a geodesic path 
on the surface and continually sheds energy into the shadow region. As 
the energy flows around the surface, it is continuously diffracted along 
the geodesic tangent toward the field point such that the significant 
effect of the surface is associated with a region around the source. In 
fact, for an ellipsoid, the significant portion of the surface, which is 
associated with the dominant energy, may look as shown in Figure 3.1. 
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ANTENNA 



REGION OF 
SIGNIFICANT 
GEODESIC PATHS 


Figure 3.1. The region of significant energy flow from an antenna 
mounted on an ellipsoid. 
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Note that this region can be specified by following the various 
geodesic paths until the radiation level along a given path becomes 
insignificant, i.e., more than 40 dB below the source magnitude. With 
this in mind, it is clear that one could represent the ellipsoid by a 
structure which simulates the elliptic cross-section completely; 
however, the profile could be approximated by a simpler shape since the 
significant energy region does not cover a large portion of the profile 
shape. An elliptic cone model is employed here to simulate the 
ellipsoid which in turn can be used to model a fuselage. This elliptic 
cone model is illustrated in Figure 3.2(a) for a source located near one 
end of the ellipsoid. Note that if the source is placed at the 
mid-section of the ellipsoid, the elliptic cone actually becomes a right 
elliptic cylinder. 

Since the elliptic cone is a developed surface, one can unfold the 
elliptic cone such that a planar structure results. Then the 
geodesics associated with the elliptic cone are straight lines on this 
planar structure. In order to allow for a geodesic solution between the 
simplicity of the elliptic cone and the rigor of the ellipsoid, one can 
perturb the elliptic cone by bending it along its generator as 
illustrated in Figure 3.2(b). In that a perturbation technique is 
employed, the geodesic paths for the elliptic cone are simply modified 
such that the solution for the ellipsoid is basically straight-forward. 
It is obvious that one cannot use this perturbation technique if 
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SIGNIFICANT 

SURFACE 

REGION 



(a) Elliptic cone simulation. 


SIGNIFICANT 

SURFACE 



(b) Elliptic cone perturbation model. 
Figure 3.2. Elliptic cone perturbation. 
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significant energy propagates far away from the source. However, as 
mentioned previously, the energy which propagates great distances along 
the ellipsoid surface becomes insignificant in magnitude such that one 
need not solve for the true geodesic paths outside the significant 
region shown in Figure 3.2(b). The simplicity of these perturbed 
geodesic paths allows one to very efficiently determine the significant 
ray paths on the ellipsoid. Thus, the efficient numerical technique, 
which uses an elliptic cylinder and elliptic cone perturbation method to 
find the creeping wave geodesic paths on an ellipsoid £13,14,22,23], 
will be discussed and verified by comparing it with the exact solution 
in Section III-C. 

It is obvious that if one wishes to simulate an aircraft, one must 
allow vertical or horizontal stabilizers, etc. to attach to the 
fuselage. The algorithm to find the curved junction edge resulting from 
attaching these plates to the ellipsoid is discussed in Section III-D. 

B. COMPOSITE ELLIPSOID MODEL 

The composite ellipsoid is used to model the general aircraft 
fuselage in this study. The composite ellipspoid is constructed from 
two ellipsoid sections positioned back to back and connected together 
such that its surface is continuous and smooth at the cross-section of 
the source location as shown in Figure 3.3. The composite ellipsoid 
semi -major/mi nor axes are defined by a 0 » b 0 , c o anc * ^o as s ^ own ™ 
Figure 3.3(a) and the source location is defined by z s and <|> s . It is 
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(a) Composite ellipsoid model. 



(b) Composite ellipsoid geometry. 
Figure 3.3. Composite ellipsoid. 
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assumed here that both the right and left ellipsoid coordinate systems 
are coincident. Then, the right side ellipsoid semi -major/mi nor axes 
and the source location are defined as 

( a f >hf ,Cf ,V gs ,V rs ) - ( a o»* 5 o* c o*^es , 'Vs^ (3.1) 

where 

.1 z s a o sin< t>s 

V es = sin c 0 and v r s " b 0 cos* s • (3.2) 

Now, one needs to define the parameters for the left side ellipsoid in 
order to complete the composition. Since the surface of the composite 
ellipsoid is continuous and smooth at the connection junction, the 
geodesic tangents of both the right and left ellipsoids should be 
continuous at the source location. For simplicity, let us consider the 
connection in the the x-z plane first. Thus, 


t e (a f ,b f ,c f ,V e$ ) = t e (a f ,b f ,c f ,V e$ ) 


(3.3) 


with 


t e ( a f. b f. c f* v es) 


-x af sinV es + z Cf cosV es 


(4 


2 22 , 1/2 

*'»»«, «s\ s ] 


and 


(3.4) 
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(3.5) 


A I I A I | 

*1,11 i i -x af sinV es + z c f cosV es 

t e (a f ,b f ,c f ,V es ) = , 2 * , , 2 2 , 7X72 

[a f sinV es +c f cos V e$ ] 


From Equations (3.3) - (3.5) one obtains 


— tan V__ 



es 


M f \ 

\J.OI 


with 


c f a d o + z sh (3.7) 

where z s ^ is the distance between the right and left ellipsoid 
coordinate origins as shown in Figure 3.3(b). Since the rim of the 
elliptic cross section at the source location is common to both the 
right and left ellipsoids, one can obtain the following relationships; 

a s ■ a f cosV es * a f cosV es (3.8) 


and 


b s * b f cos V e$ - b f cosVg S 


(3.9) 
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From Figure 3.3(b), one also finds that 


c f sinV es = 2 $ - 2$h . 


(3.10) 


Now, from Equations (3.7) and (3.10) one obtains 


z sh^ s1nV es + U = z s - d o sinV es 


(3.11) 


i 


* ® • * 


i 

2 S - d 0 si nV es 

z sh = 1 + sinV* 

es 


(3.12) 


After manipulation of Equations (3.3) - (3.12), the parameters needed to 
define the left ellipsoid are given by 


I 


V 


es 


. -1 
sin 


|- c o cos v es 

L tanV es {d Q + z $ ) * 



(3.13) 


' a o cos v es 

a f ” cos V’ 
es 


• b o cos v es 

b f * cos V- 
es 


i 

and Cf was already given in Equation (3.7). 


(3.14) 


(3.15) 
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The smoothness through the whole elliptic rim at the connecting 
junction can be easily checked by applying Equations (3.6), (3.7), 

(3.14) and (3.15) to the x, y, and z components of t g and n, which will 

be given in Equations (3.24) and (3.26), respectively, for an arbitrary 
value of V r . It is also noted that the two surfaces join together at 
the source location so creeping waves that propagate to the right side 
from the source will follow the geodesics on the right ellipsoid surface 
and those that propagate to the left side from the source will follow 
the geodesics on the left ellipsoid surface. 

C. NUMERICAL TECHNIQUES FOR THE GEODESIC CALCULATION OF AN ELLIPSOID 

It is seen that, for an antenna mounted on an ellipsoid, the 
geodesic paths associated with the UTD solution in the shadow region are 
extremely complex. An elaborate method employing calculus of variations 
to calculate the geodesic paths, which is very inefficient, is presented 
in Appendix A. An efficient numerical approach is examined in this 
section with the ellipsoid simulated by an elliptic cone or elliptic 
cylinder model. Since the elliptic cone and elliptic cylinder are 
developed surfaces, geodesics are efficiently obtained. 
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1. Surface Geodesics 


The geodesics on the elliptic cylinder, elliptic cone and ellipsoid 
are presented in this section. 

(a) Elliptic Cylinder Case: 

The elliptic cylinder geometry used for this study is shown in 
Figure 3.4(a). Since the elliptic cylinder is a developed surface, the 
geodesic path TH} is a straight line on the unfolded planar surface as 
shown in Figure 3.4(b). Thus, the geodesic unit tangent t is given by 

t « V r COSY + 2 SinY ( 3 . 16) 


where 


a * 

x(-a sinV r ) + yb cosV r 

, 2.2 2 2 1/2 

(a*sin^V r + b cos v p ) 


(3.17) 


It is noticed that V r and z are the two principal directions on the 
elliptic cylinder surface. For a given geodesic Q ' Q , one can see that y 
is a constant along the geodesic path. 
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SOURCE 


A 
!• 



(a) TRUE ELLIPTIC CYLINDER 



(b) UNFOLDED PLANAR SURFACE 


Figure 3.4. Geodesic path on a developed elliptic cylinder. 


(b) Elliptic Cone Case: 


Consider a ray which propagates along a geodesic path THT on the 
elliptic cone surface as shown in Figure 3.5(a). It is a straight 
line on the unfolded planar surface as shown In Figure 3.5(b). It is 

* A 

noticed that t g and tj are the two principal directions on the surface. 

A 

The goedesic unit tangent t is, then, represented by 



cosB + 


sin 8 


where 

t g = -x g sin 6 + z cos 6 


(3.18) 


(3.19) 


A A 

xa< cosV r + yb* sinV r 
< /a'^cos 2 V r + b'sin 2 V r 


(3.20) 


and 6 is the half cone angle as shown in Figure 3.5(a). Note that 8 

* A 

denotes the angle between t^ and the geodesic unit tangent t and is no 
longer a constant along the geodesic path Q'Q. 


(c) Ellipsoid Case: 

Recall that the ellipsoid surface is defined by 


R(V e , v r ) 


a cosV e cosV r x + b cosV g sinV^ + c sinV g z 


with 


(3.21) 


SOURCE 



(a) TRUE ELLIPTIC CONE 



(b) UNFOLDED PLANAR SURFACE 
Figure 3.5. Geodesic path on a developed elliptic cone. 
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tan V. 


cot 0 


- sin 2 4 
_~b2~ 


2 1/2 
COS £ «T| 


a 2 J 


(3.22) 


and 


tan V„ « 
p 


a si n 
b cost 


(3.23) 


As shown in Figures 3.6(a) and (b), the unit vectors in the principal 
directions can be represented by 


3R 



-xa sinV e cosV r - yb sinV e sinV r + zc cosV e ^ 

[a 2 sin 2 V cos 2 V +b 2 sin 2 V sin 2 V +c 2 cos 2 V 1 1/2 
e r e r e 1 


and 


t 


1 = t e x n 

-xa sinV r (b 2 sin 2 V e + c 2 cos 2 V e ) + yb co sV r (a 2 sin 2 V e 

[a 2 b 2 sin 2 V +c 2 cos 2 V (a 2 sin 2 V +b 2 cos 2 V )] 1/2 
i e e r r 1 


+ c 2 cos 2 V e ) 


+ zc(b 2 -a 2 ) sinV r cosV r sinV e cosV e 
[c 2 cos 2 V g +sin 2 V e (a 2 cos 2 V r +b 2 sin 2 V r ) ] 1/2 (3.25) 
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(a) Geodesic path on an ellipsoid. 



*• VO 2 cos 2 V, + b 2 sin 2 V, 


(b) Elliptic profile. 

Figure 3.6. Definition of the unit vectors on an ellipsoid. 
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with 


t r x t e 

n = — — 

|t r x t e | 


and 


xbc cosV e 
[a 2 b 2 sin 2 V 

e 


cosV r + yac cosV e sinV r + zab sinV e 
+ c 2 cos 2 V e (a 2 sin 2 V r + b 2 cos 2 V p ) } 1 ' 2 




9R 




-xa cosV e sinV r + yb cosV e cosV r 
cos 2 V g sin 2 V r + b 2 cos 2 V e cos 2 V p ] 


1/2 


(3.26) 


(3.27) 


If S denotes the angle between tj and the geodesic unit tangent t, then 
* 

t is represented by 


t 3 tj cos$ + t fi sine ^3 2 

which is identical to the form used for the elliptic cylinder and 
elliptic cone geodesics. This suggests that one might be able to 
develop a perturbation solution which gives a simplified form for e on 
an ellipsoid using the elliptic cylinder or elliptic cone expressions 
for 6. 
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2. Elliptic Cylinder Perturbation 


In order to solve for the geodesics on an ellipsoid, the elliptic 
cylinder perturbation technique is used when the source is located at or 
near e s = 90°. Recalling that y is a constant along a given geodesic 

path on an elliptic cylinder, one obtains a geodesic equation which is 
given by 


tan y 


Se 


(3.29) 


The elliptic cross-section in the x e - z plane is described by 


x_ = 


/ 


2 2 
a^cos^V. 


b^sin^V. 


cosV. 


(3.30) 


and 

z = c sinV e 


The arc length in the t e direction is given by 


(3.31) 


V e I 3x^ 2 3z r 

s e = £ /(«;) ♦ dY e 



(3.32) 
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The elliptic cross-section in the x-y plane is described by 
x = a cosV g$ cosV p , and 

y ■ b cos V gs sinV p 

The arc length in the t r direction is given by 


S r = 


f 

V r 


dx d dy 

tin) + tin) 


dV. 


» p o o i o 0*1/2 * 

= / cosV es {a^sin^V r + b^cos^} dV p 

Vr _ 


Thus, the geodesic equation becomes 


tan y * 


Ve 9? 9' 9 9'^^ ' 

/ [(a cos V p + b sin^V p ) sin^V g + c^cos^V e l dV £ 


V r? ? 0 1 0 O' 0 • 1/2 , 

/ [a^cos^V es sin^V r + b‘cos‘v es cos\] dV p 


Assuming the diffraction point is located at Q (a cosV g cosV 
b cosVg sinV r , c sinV fi ) and the field point at P (R t sine^ cos<j> t , 
R t sin0 t s i n <f> t , R t cos9 t ), then at the diffraction point Q, the 


( 3 . 33 ) 

( 3 . 34 ) 


( 3 . 35 ) 


( 3 . 36 ) 
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radiation direction should coincide with the geodesic tangent t 

given in Equation (3.28). Thus, the radiation direction is given by 

t = xt x + + zt z 

= tj COST ♦ t e Sin, , ( 3 . 37) 


where 


t 


x 


a 

sine t cos <j» t - j- cosV e cosV r 


D 


(3.38) 


y 


b 

sine t si n <t> t - cosV g cosV p 


D 


(3.39) 


and 


t 


z 


c 

cos sinV e 


D 


(3.40) 


Note that 




b cosV p cosy(a 2 sin 2 V e + c 2 cos 2 V g ) 

[a 2 b 2 sin 2 V e + c 2 cos 2 V e (a 2 sin 2 V p + b 2 cos 2 V r ) ] 1/2 

1 

[c 2 cos 2 V e + sin 2 V e (a 2 cos 2 V r + b 2 sin 2 V p ) ] l/2 
b sinV & sinV r siny 

[c 2 cos 2 V e + sin 2 V e (a 2 cos 2 V p + b 2 sin 2 V p ) ] 1/2 
b 

sine t si n 4> t - cosV g sinV p 

D » (3.43) 


c(b 2 -a 2 )sinV p cosV p sinV fi cosV g cosy 

[a 2 b 2 sin 2 V e + c 2 cos 2 V e (a 2 sin 2 V p + b 2 cos 2 V p ) ] W2 

1 

[c 2 cos 2 V g + sin 2 V e (a 2 cos 2 V p + b 2 sin 2 V p ) f' 2 

c cosV e siny 

[c 2 cos 2 V g + sin 2 V e (a 2 cos 2 V p + b 2 sin 2 V p ) ] W2 
c °s0 t -J- sinV g 


D 


(3.44) 



Now, ft„b cosV„ + t,a sinV„1 c cosV a + t, ab sinV. yields 

L x r y r J e z e 


H( ®t» *t, V e , V r , R t ) 

= absinV e cose^ + csin8 t cosV g (asi n <t> t sinV p + bcos<f> t cosV p 
abc 

- w~ = 0 • 

Next, from Equations (3.42) and (3.43) one obtains 

t x (-b sin V p ) + t y (a cos V p ): 

ab cosy{sin^V e (b 2 sin 2 V r + a 2 cos 2 V r ) + c 2 cos 2 V e } 1/2 
{a 2 b z sin 2 Vg + c z cos 2 V e (a^si n i; V r + b 2 cos 2 V p ) } 1/2 

= q" {acosV p sin0 t si n - bsinV p sin0 t cos$ t } . 
Accordingly, 

1/2 

cosy = [a 2 b 2 sin 2 Vg + c 2 cos 2 V e (a 2 sin 2 V r +b 2 cos 2 V r ) ] 

sin6t(acosV r sin<^-bsinV r cos^) 
Dab{sin 2 V e (b 2 sin z V r +a z cos 2 V r ) + c 2 cos z V g } 1/2 


(3.45) 


(3.46) 


(3.47) 
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Substituting Equation (3.47) into Equation (3.44), one obtains 


[c(b 2 -a 2 )SinV p 'cosV p sinV e cosV e sin9 t 

• (a cosV p sin$ t - b sinV p cos$ t ) 

Se 

+ c cosV fi 3j: sine t (a cosV p sin<j> t - b sinV p cos$ t ) 

• {a 2 b 2 sin V e + c 2 cos 2 V e (a 2 sin 2 V r + b 2 cos 2 V r )} W *] 

ab[c 2 cos 2 V e + sin 2 V e (a 2 cos 2 V p + b 2 sin 2 V p ) ] 

c 

* cos 9 t - sin V e 


(3.48) 


Thus, one finds that 

G ( e t ,<f> t ,V e ,V p ,R t ) = S p c(b 2 -a 2 ) sinV p cosV p sinV g cosV g sin9 t 

• (a cosV p sin$ t - b sinV p cos<fr t ) 

+ S g c cosV e [a 2 b 2 sin 2 V e + c 2 cos 2 V e (a 2 sin 2 V p + b 2 cos 2 V p ) ] 1/2 

• sin9 t (a cosV p si n 4> t - b sinV p cos<j> t ) 

- S p ab cos 9^ [c 2 cos 2 V g + sin 2 V e (a 2 cos 2 V p + b 2 sin 2 V p )] 

s 

+ abc sinV e [c 2 cos 2 V g + sin 2 V e (a 2 cos 2 V p + b 2 sin 2 V p )"| 


* 0 


(3.49) 
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Note that S e and S r were given in Equations (3.32) and (3.35), 
respectively. 

Provided that one has obtained a diffraction point (V e , V r ) for a 
receiver location (Rt» » a numerical technique can now be developed 

from Equations (3.45) and (3.49) to solve for (V e + AV e , V r + AV r ) 
associated with a new receiver location (R^ + AR t , 9t + A9 t> + 
Assuming that the i th set of (R t , e t , V e , V r ) is first known to 
satisfy Hj = G-j ■ 0, or at least approximately so, the next set 
(R t + AR t , e t + A9 t , ^ v e + AV e , V r + AV r ) is obtained by 

enforcing H-j+j = G-j+j = 0, such that 


« 1+1 ■ H i + H R t 4R t * H e t 49 t + \‘*t * 


H v 4V e t H v lV p 
e r 


= 0 


and 


(3.50) 


6 1+1 • G, + G„ t 4R t * ♦ G^4» t ♦ 


t G v iV + Gy 4» r 
e r 


* 0 


(3.51) 
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In matrix form, these equations are satisfied 



V 


AV e 


~- H 1 - - H^M t - 

\ 

GVr _ 


AV r 


-S, - 


Note that the partial derivatives are given by the following: 


Hy fi = ab cosV e cos9 t - c sinV e sin9 t (asin<j> t sinV r + bcos<f> t cosV r 
Hy r = c cosV e sin9 t (asin<}> t cosV r - bcos<}> t sinV r ) 

H 0 = -ab sinV fi sin9 t + c cosV e cos9 t (asin<f> t sinV r + bcos^cosV 


- c cosV fi sin9 t (a cos<|i t sinV r - b sin$ t cosV r ) 
t 

^R t = abc/R^ 


( 3 . 52 ) 

( 3 . 53 ) 

( 3 . 54 ) 

) 

( 3 . 55 ) 

( 3 . 56 ) 

( 3 . 57 ) 
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G v = S p c(b 2 -a 2 ) sinV r cosV r (cos 2 V e -sin 2 V e ) 

• sin0 t (a cosV p sin<t» t - b sinV p cos<j> t ) 
dS 

+ 3 ^ c cosV e [a 2 b 2 sin 2 V e + c 2 cos 2 Vg(a 2 sin 2 V p + b 2 cos 2 V p )] 

• sine t (a cosV p sin<(> t - b sinV p cos<|> t ) 

+ S g c sin0 t (a cosV p sin<j> t - b sinV p cos$ t )sinVg 

a 2 b 2 (cos 2 V -sin 2 V ) - 2c 2 cos 2 V (a 2 sin 2 V + b 2 cos 2 V ) 
e e e r r 

[a 2 b 2 sin 2 V + c 2 cos 2 V (a 2 sin 2 V + b 2 cos 2 V )] 1/2 

L g p p/ j 

+ 2 S p abcos0£ cosV e sinV g (c 2 -a 2 cos 2 V p - b 2 sin 2 V p ) 

S p 

+ ^ abc cosV e [c 2 cos 2 V e + sin 2 V e (a 2 cos 2 V p + b 2 sin 2 V p ) 

+ 2sin 2 Vg(a 2 cos 2 V p + b 2 sin 2 V p - c 2 )] 


where 


^ = {c 2 cos 2 V e + (a 2 cos 2 V p + b 2 sin 2 V p )sin 2 V e } 1/2 
e 


1/2 


(3.58) 


(3.59) 
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and 


, u,s r 9 9 

G V p * 1 "3V7 sinV p cosV r + S p (cos z V p - sin z V p ) } 

• c(b 2 -a 2 ) sinV e cosV e 'sin 9 t (a cosV p sin<f> t - b sinV p cos* t ) 
- S p c(b 2 -a 2 ) sinV p cosV p sinV e cosV fi sin 0 t 

• (a sinV p sin<t» t + b cosV p cos<f> t ) + 



• sine t (a cosV p s i n <t> t - b sinV p cos^) 

S e c 3 cos 3 V e (a 2 -b 2 )sinV p cosV p sin 9 t (acosV p sin$ t -bsinV p cos4> t ) 
[a 2 b 2 sin 2 V g + c 2 cos 2 V e (a 2 sin 2 V p + b 2 cos 2 V p ) ] 1/2 

" s e c cosV e [a 2 b 2 sin 2 V g +c 2 cos 2 V e (a 2 sin 2 V p + b 2 cos 2 V p ) ] 1 ' 2 
• sin9 t (a sinV p sin<j> t + b cosV p cos<|> t ) 
dS 

- ab cos9 t [c 2 cos 2 V e + sin 2 V e (a 2 cos 2 V p + b 2 sin 2 V p )] 

+ 2S p ab(a 2 -b 2 ) cos9 t sin 2 V g cosV p sinV p 

dSp abc + 99 999 99 

+ sinVgfc^cos^Vg + s1n^V e (a^cos^V p + ITs1rrV p )] 

2S 

+ abc sin 3 V e (b 2 -a 2 ) cosV p sinV p (3,60) 
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where 


dS 


UJ r , ? o ? ? ,1/2 

* cosV es ( aZsin V + b - cos v r ^ 


(3.61) 


and 


d^e Ve 

HV7 = / 

0 


(b 2 -a 2 ) sinV cosV sin 2 v‘ 
r r e 


[c 2 cos 2 V a + (a 2 cos 2 V_ + bSin^Vjsin^V-V 
l e r r' e 1 


2.,.2„ ,.,.2„V /2 dV ( 


(3.62) 


and 


G ft = S„ c(b 2 -a 2 ) sinV cosV sinV cosV cose* 

0^ r v r r e e t 

• (a cosV p s i n <)> t - b sinV p cos^) 

+ S g c cosV e [a 2 b 2 sin 2 V e + c 2 cos 2 V e (a 2 si n 2 V p + b 2 cos 2 V p )1 1/2 


• cose t (a cosV p s i n - b sinV p cos 4> t ) 

+ S p absin0 t [c 2 cos 2 V e + sin 2 V e (a 2 cos 2 V p + b 2 sin 2 V p )] 


(3.63) 


= S p c(b 2 -a 2 ) sinV p cosV p sinV e cosV e sin0 t 

• (a cosV p coS(j» t - b sinV p sin$ t ) 

+ S e c cosV e [a 2 b 2 sin 2 V e + c 2 cos 2 V e (a 2 sin 2 V p + b 2 cos 2 V p )] l/2 

• sine t (a cosV p cos^ + b sinV p sin<f> t ) (3.64) 
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and 


5 r 

G r = - rT abc sinV e [c 2 co$ 2 V e + sin 2 V e (a 2 cos 2 V r + b 2 sin 2 V )] . 
t t ^ 

(3.65) 


It is seen that one can solve for (AV e> AV p ), for a known (AR t ,A0£,A<j> t ) , 

IlCinn 4 am f O CO\ T* £ i.i • * • i j . i •• . • . 

u M uauiun iu ueteriii i ne une initial avnracnon point 

(V e ,V p ) for a given receiver location (R t ,e t ,$ t ), one can always assume 
a diffraction point at the source (0 »V ps ) with the radiation direction 
(9f»4y = ±ir/2) with respect to the source coordinate system (t N ,t p ,t e ) 
and gradually add increments (AR t ,A0 t ,A<f> t ) until the desired receiver 
location ( R t » is reached as depicted in Figure 3.7. 

In this process one can construct a cone where the rim of the cone 
is traced out by the receiver trajectory with the tip of the cone at the 
source point Q' and the cone axis aligned with t g . The half cone angle 
is given by 


9^ * tan 


fi 


^N ,p osl 2 + IV p osl 2 


V p os 


(3.66) 


with 


P 0S * P 0 " P S (3.67) 
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/ PROJECTION TO \ 
\ t N - t r PLANE / 


Figure 3.7. Illustration of the diffraction point finding for a 
given receiver location. 
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and 4y(N) is given by 


4y(N) = tan"* 


VPgs 

V P OS 


(3.68) 


Note that P(l) denotes the position vector of the assumed observation 
point tangential to the source, i.e., (9^.,^ = ir/2 ) with respect 

AAA 

to (t^,t r ,t e ) and P(N) denotes the position vector of the actual 
observation point tangential to the diffraction point Q, i.e., 

(9 f ,<fr f + NA$ f ) with respect to (t N ,t p ,t e ) or (9 t ,4> t ) with respect to 

AAA 

(x.y.z). It is observed that there exists a one-to-one correspondence 
between the points (from P(l) to P(N)) on the rim of the cone and the 
points on the ellipsoid surface. 

After the initial diffraction point is identified by (V e ,V r );y, 
one knows the geodesic path as defined by the geodesic equation which 
was given by 

S e 

tan ^ = 37 (3.69) 

and y is a constant along a given geodesic path on the perturbed 
elliptic cylinder. Such a numerical approach is illustrated in Figure 
3.8. One need not trace out the complete geodesic path from the source 
location to the diffraction point for each new radiation direction. As 
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X 



( 6 f +A0 1 i ) 


Figure 3.8. Elliptic cylinder perturbation. 
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shown in Figure 3.8, the diffraction point (V e + AV e , V r + AV r ) for the 
next receiver location is determined from (V e ,V r ), using Equation 
(3.52), if (AR t ,A0 t ,A(^) is small which is the case when a complete 
radiation pattern is computed. 

After the geodesic path is determined, various other parameters 
associated with actual field calculations must be found, the Fock 
Darameter E of Fnuat.inn ic rrivon hv 

* * * ' ■ ' " » X— ---/ ' - 3 ’ « - —j 

Q 1 kp g 1/3 

C = ^ "^g ( 2 ) d * * (3.70) 


where the integral is evaluated along the geodesic path. Note that t is 
the arc length along the geodesic and given by either of the following 
equations; 


l = 


siny 


or l = 


cosy 


(3.71) 


Accordingly, 


d* 


siny 


or 


/ 


(a 2 cos 2 V r + b 2 sin 2 V p ) sin 2 V g + c 2 cos 2 V g dV £ 


(3.72) 


d* 


J_ / 

cosy l 


a 2 cos 2 V es sin 2 V r + b 2 cos 2 V es cos 2 V p dV p ( 3 . 73 ) 
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and the integration can be rewritten as 
1 V r 1 k 0g 1/3 . 

5 ■ Si 7 f ■'J (t) 

Vr s 


•/ a 2 cos 2 V es sin 2 V r + b 2 cos 2 V es cos 2 V p dV p (3.74) 


or 


5 = 


siny 


r 

0 




c 


ko 1/3 

p £) 


• /( a 2 cos 2 V + b 2 sin 2 VJ sinV + c 2 cos V dv‘ 
j r r' e e e 


(3.75) 


2 2 

where = l/(K^cos y + KgSin y) and K^, are two principal curvatures 
given in Equations (2.79) and (2.80), respectively. 


Next, the ray divergence factor )/d<|>(Q) is defined as 

the change in the separation of adjacent surface rays as shown in 
Figure 3.9. Since the ellipsoid simulating the aircraft fuselage will 
be long and slender, it is assumed that the ray divergence factor is 
unity in the analysis. 

This completes the elliptic cylinder perturbtion solution for an 


antenna mounted on the mid-section of the ellipsoid. 
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3. Elliptic Cone Perturbation 

When the source is not located at or near the mid-section (0 S =9O°), 
the elliptic cone perturbation method is used. The geodesic path on an 
ellipsoid using the elliptic cone perturbation method is shown in Figure 
3.10(a) and the associated unfolded surface is shown in Figure 3.10(b). 

A A 

If y and 5 denote the angle between t and t 1 at Q' and Q, respectively, 
it was shown in Section III-C.l that they are not the same as in the 
elliptic cylinder case. In fact, 

8 - Y - a. (3.76) 

The calculation of a will be discussed in detail in this section. 

Recall that the ellipsoid surface is defined by 

R(V e ,V p ) = a cosV e cosV p x + b cosV g sinV^ + c sinV e z (3.77) 

-► 

and the tangent vector t e is defined by 

-► 3R a a ^ 

t e * s -xa sinV e cosV p - yb sinV e sinV p + zc cos V g . (3.78) 

Therefore, one can show that 
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SOURCE 


DIFFRACTION 

POINT 


Wfi 



(a). Geodesic path on the perturbed elliptic cone. 



(b). Geodesic path on a developed elliptic cone. 

Figure 3.10. Illustration of the elliptic cone perturbation method. 
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(3.79) 


x-a cosV es cosV rs y-b cosV es sinV rs z-c sinV es 

-a sinV es cosV rs * -b sinV es sinV rs * c cosV es 


From Figure 3.11, z * c esc V es for x = y = 0; thus one obtains 


[a 2 cos 2 V cos 2 V ps + b 2 cos 2 V fiS sin 2 V ps 


+ (c sinV es - c esc V gs ) ] 


1/2 


• [.* cos 2 V rs ♦ b 2 s,n 2 V rs ♦ z 2 cot\ es 1 


1/2 


where 


z $ * c sinV e$ 


a = a cosV 

a$ a 


b = b cosV 
$ es 


From Figures 3.12(a) and (b). 


% * x a $ cosV r + y b $ sinV r 


and 


* 

at 


» - x a $ sinV p + y b $ cosV p 


(3.80) 


(3.81) 


(3.82) 
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Figure 3.11. Elliptic profile. 
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From the previous calculation, one can define a new vector U as 
follows: 


-► •¥ 

* at aR 

u * CaV^) x CaV^l 

A A A 

(- x a sinV gs cosV r - y b sinV gs sinV r + z c cosV e$ ) 

A A 

• (- X a s sinV r + y b,.cosV p ) 

[a 2 sin 2 V e$ cos 2 V p + b 2 sin 2 V e$ sin 2 V p + c 2 cos 2 V gs 1 1/2 « 

O O O O 1/2 

• [a^sin V p + b*cos V p ] 



( 3 . 83 ) 
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The magnitude of vector U becomes 


|U| 


-C" 

1L. 

3V r 


9R 


sine 0 = s1ne n 


From Figure 3.12(b), one obtains the following equation; 


da = 


ds r |U| 


Therefore, a is given by 


a = 


V 

i 


r 


V 


rs 


cosV es [a 2 b 2 sin 2 V es + c 2 cos 2 V es • (a 2 sin 2 V r + b 2 cos 2 V r )] dV 
r[c^cos^V es + sin^V e$ (a 2 cos^V^ + b^sin^V^.) ] 1/2 ~ 


where 


rs 


Lafbf + z 2 cos 4 V es (a 2 sin 2 V^. + b 2 cos 2 v|.) ] 1/z 
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a^cos^V^ + b^sin^Vj. 


+ z^cot^V 


es 


dV. 


. 2 2 2 2 2 4 , 1/2 

r * [a‘cos^V r + b^sin ^V p + z*cot*V es ] 


(3.85) 


(3.86) 


(3.87) 
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To obtain a diffraction point (V e ,V p ) for a radiation direction (e t ,4> t ), 
one follows the same idea as that of the elliptic cylinder perturbation 
solution. 

At the diffraction point Q, the radiation direction ( © t , ) should 

A 

coincide with the geodesic tangent t. Thus, 

A A A A 

t = X t x + y t y + z t z 

* t x cos e + t e sin 0 t (3.88) 

Note that t , t , t , t , and t. were already given in the Equations 

A Jr X C X 

( 3 . 38 ) — ( 3 . 40 ) and (3.24) and (3.25), respectively. Equating the x-, y-, 
and z-components , respectively, one obtains 

-a sinV p cos$(b 2 sin 2 V e + c 2 cos 2 V e ) 

[a 2 b 2 sin 2 V + c 2 cos 2 V (a^sin^V + b^cos 2 V )] 1/2 
L e e r r 

i 

[c 2 cos 2 V g + sin 2 V g (a 2 cos 2 V r + b 2 sin 2 V p ) ] 1/2 

asinV^ cosV„ sins 
e r 

• 9 9 9 9 9 9 9 .1/2 

[c^cos^V e + sin V fi (a^cos v p + b^sin^V p ) ] 
a 

sine t cos$ t - cosV g cosV p 

= D » (3.89) 
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b cosV p cosB(a 2 sin 2 V g + c 2 cos 2 V g ) 


[a 2 b 2 sin 2 V g + c 2 cos 2 V' g (a 2 sin 2 V r + b 2 cos 2 V r ) ] 1/2 

1 

[c 2 cos 2 V g + sin 2 V g (a 2 cos 2 V r + b 2 sin 2 V p ) ] 1/2 


b sinV. sinV_ sinB 
e r 


[c 2 cos 2 V g + sin 2 V g (a 2 cos 2 V p + b 2 si n 2 V p ) ] 1 


and 


sine t cos$ t - cosV e sinV p 


( 3 . 90 ) 


o o 

c(b - a ) sinV cosV„ cosV cos$ 
v ' r r e 




[a 2 b 2 sin 2 V fi + c 2 cos 2 V e (a 2 sin 2 V r + b 2 cos 2 V p ) ] 1/2 

1 

[c 2 cos 2 V g + sin 2 V e (a 2 cos 2 V p + b 2 si n 2 V p ) ] 1/2 


c cosV e sinB 


[c^cos^Vg + sin*V e (a cos V p + b^sin^V p ) ] 


cose* - s1n\L 
t R t e 


( 3 . 91 ) 
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From Figure 3.10(b), one can show that the perturbed geodesic path can 
be expressed as follows: 


r g cos(y-a) * r $ cosy 


where 


r. 

s 


fa 2 cos 2 V 
- s rs 


+ b?sin 2 V__ + z?cot 4 V^, 

* r a a ci 


1/2 


(3.92) 


(3.93) 


r = (a 2 cos 2 V 


2 2 2 4 n 1/2 

r + b^sin V r * z‘cot\ s ) - S ( 


(3.94) 


and 



2 2 ' 

[c cos V + 


2 2 
(a cos V 


+ b 2 sin 2 V p )sin 



1/2 , 
d\L 


(3.95) 


Now, [t x (-b sinV r ) + t y a cosV p ] yields 

ab cos(y-o) [c 2 cos 2 V e + sin 2 V e (a 2 cos 2 V r + b 2 sin 2 V r )] 
[a 2 b 2 sin 2 V e + c 2 cos 2 V e (a 2 sin 2 V p + b 2 cos 2 V p ) ] 1/2 


sin e t (acosV r sin^ - bsinV r cos^) 
D 


(3.96) 
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Next, [(t x b cosV p + t a sinV p ) c cosV g + t z a b sinV e ] yields 


absinV e co$9 t + csine t cosV e (asin^sinVp + bcos$ t cosV r ) 

_ - 

abc 

- DR7 = o . (3.97) 

Three functions can, then, be constructed from Equations (3.92), (3.96), 
and (3.97) as follows: 

F(V e ,V p ,Y) = r e cos(y-a) - r $ cosr = 0 t (3.98) 

G(R t . 8 t> * t . »,. » r . 1 ) 

= DabcostY-aJCc^cos^Vg + sin^Vgta^cos^Vf + b2sin^v r )] ^ 

- sin0 t ( a s i n 4> t cosV p - bcos<t> t sinV p ) 

• [a 2 b 2 sin 2 V e + c 2 cos 2 V e (a 2 sin 2 V r + b 2 cos 2 V r )] 1/2 
= o 0.99) 

and 

H(R t , 6 t , * t . V e , V r ) 

= absinV e cos9 t + csin0 t cosV g (asin^ sinV p + bcos<t> t cosV p ) 
abc 

- » 0 . (3.100) 
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Provided that one has obtained a diffraction point (V g , V p ) for a 

receiver location (R t , © t , 4> t ), a numerical technique can now be 

developed from Equations (3.90) (3.99) , and (3.100) to solve for 

(V + AV . V + AV ) associated with a new receiver location (R t + AR t , 
v e e r r c L 

9^ + A0 t> + A<j> t ). Assuming that the i th set of (R t , 0^, $ t , V g , V p ) 
is first known to satisfy F 1 = ^ = G i = 0, or at least approximately 
so, the next set (R t + AR t , 0 t + A0 t , + A$ t , V 0 + AV e , V p + AV p ) is 
obtained by enforcing = ^i+i = suc * 1 

F i+1 * F i + F V AV e + F V AV r + F y Ay " 0 » (3.101) 

e r 


G i+1 = G i + G v AV e + G V AV r + G y Ay 
e r 

+ G Q A0 t + G^ A<(> t + G r AR t = 0 
t t t 


(3.102) 


and 


H 1n * H ( + H v 4V e ♦ H y 4V p ♦ H, H 
6 r t 


+ A«j» t + H r AR t s 0 
t t 


(3.103) 
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In matrix form, these equations are satisfied by 


— — 


— — 



F v F v r f y 


AV e 


■- F i 

e r 





G V Gy G 

e r T 


AV r 

- 

' G i - G e t 49 t ‘ G * t 4 *t - G R t 4R t 

Hy Hy 0 
e r 


ay 


- H 1 - - H^ t - 


(3.104) 

Note that the partial derivations are given by the following: 

Fy = -[c 2 cos 2 V e +(a 2 cos 2 V r +b 2 sin 2 V r )sin 2 V e ] 1/2 cos(y-a) , (3.105) 


F v = (b 2 -a 2 ) sinV p cosV r 
r 


cos V 


es 


(a^cos^V r +b^si n^V r +z^cot 4 V e$ ) 1 /2 


V e sin-V e 

Jy [c^cos^Vg+U^cos'^+b^si n^V^sin^Vg ] 1 /2 



cos (y-a) 


r e si n ( y-«) [a 2 b 2 +Z 2 cot 4 V fis (a 2 si n 2 V p +b 2 cos 2 V p ) 1 1 /2 
a^cos^V p +b^si n^V p +Z^cot 4 V es 


(3.106) 
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* 


(3.107) 


F y = r $ sin y - r g sin (y-a) 


and 


Dab(a 2 cos 2 V r +b 2 sin 2 V r -c 2 )sinV e cosV e cos (Y-a) 
\rZ^r 2w TnZTZTZxi xk2 r i n 2 \/ \c i nZ\/ lT/T 


^ u J 


' L/ Jill 


'e J 


2 2 2 2 2 2 2 

sin e t (asin4> t cosV r -bcos4JtSinV r )Ca b -(a sin V r +b cos V r )c ]sinV e cosV e 

[(a^sin^V.+b^cos^V-^cos^V.+a^b^sin^V. ] 1/2 
L r r ' e 6 J 


abcOS(Y-a) 


[c 2 cos 2 V e +(a 2 cos 2 V p +b 2 sin 2 V r )sin 2 V e ] 


1/2 


c 2 -(a 2 cos 2 Vp+b 2 sin 2 V r ) 


sinV -cosv. 
e e 


+ sin9 t sinV e (R^ cos<f> t cosV p 


b c 

+ ^ sin<|) t sinV r J - cose t cos V g 


(3.108) 
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Dab co s(Y-a)(b 2 -a 2 )sinV r cosV r sin 2 V e 

G V r “ [c 2 cos 2 V e +(a 2 cos 2 V p +b 2 sin 2 V p )sin 2 V e ] 1/2 

+ Dab[c 2 cos 2 V e +(a 2 cos 2 V r +b 2 sin 2 V r )sin 2 V e ] 1/2 sin(y-a) 

[a 2 b 2 +Z 2 cot 4 V es (a 2 sin 2 V r +b 2 cos 2 V r )] 1/2 

a 2 cos 2 V r +b 2 sin 2 V r +Z 2 cot 4 V es 

sin9 t (asin4> t cosV r -bcos^.sinV r )(b 2 -a 2 )sinV r cosV r c 2 cos 2 
+ [(a 2 sin 2 V r +b 2 cos 2 V r )c 2 cos 2 V e +a 2 b 2 sin 2 V e ] i/2 


2 2 2 2 

+ sine t (asin$ t sinV r +bcos$ t cosV r ) [(a sin V p +b cos V p ) 


c 2 cos 2 V +a 2 b 2 sin 2 V] 


1/2 


abcOS(Y-a) o o 0 o ? ? ? ,1/2 

— “5 {c 2 cos v e +(a^cos v p +b^si n^V p )si n*V g ] 


(b 2 -a 2 ) 2 b__ . 

{ — sinV p cosV p cos £ V e -sin9 t cosV e (R t sin<i> t cosV p 


- ^ cos* t sinV r ) } 


O O O O O 0 2 -I 

* -Dabsin(Y-a)[c^cos^V e +(a cos V p +b sin V p )sin V e ] 


( 3 . 109 ) 


( 3 . 110 ) 
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- cos9 t (asin<fi t cosV r -bcos<|> t sinV r ) 


• [(a 2 sin 2 V r +b 2 cos 2 V r )c 2 cos 2 V e +a 2 b 2 sin 2 V e ] 1/Z 


1/2 


abcos(ir-a) 9 9 9 ? ? ? 9 , 

jj [c 2 cos 2 V e +(a z cos z V r +b z sin z V r )sin z V e ] 

c a b 

{ sin0 t sinV e -cos0 t cosV e (R t . cos^cos V f + si n 4. t si nV p J } , 

(3.111) 


• - sin© t (acosi}) t cosV r +bsin(ji t sinV r ) 

• [(a 2 sin 2 V +b 2 cos 2 V )c 2 cos 2 V +a 2 b 2 sin 2 V ] 1/2 

r r c c 

abCOS(Y-a) o o 9 9 9 9 9 , 1/2 

+ p [c 2 cos z V e +(a z cos z V r +b z sin z V r )sirrV e l 

a b 

• { sin 0 t cosV e (r^ sin 4 . t cosV r - cos<() t sinV r )} , (3.112) 


ab ,99 99 22 .? <1/2 

y cos(Y-a)[c z cos z V e +(a z cos z V r +b z sin^V r )sin^V e l 


2 

{[sin0 t cosV e (a cos<t> t cosV r +b sin^sinV,.) + c cos0 t sinV e ]/R t 


[cos 2 V e (a 2 cos 2 V r +b 2 sin 2 V r )+c 2 sin 2 V e )/R^} 


(3.113) 




H 


= ab cos V cose* - c sin\Lsine*(asin<t> t sinV„ + bcos^cosVj 

e 1 e t t r t r (3ai4) 

= c cosV e sine t (asin* t cosV r - bcos$ t sinV p ) (3.115) 

(3.116) 

» 0 

* -ab sinV e sin9 t + c cosVgCOse^asin^sinVj. + bcos<|> t cosV r ) 

(3.117) 

= c cosV e sine t (a cos$ t sinV r - b sin^ t cosV p ) (3.118) 


and 

abc 

H r = J2~ • (3.119) 

* t 

It is seen that one can solve for (AV g , AV r , Ay), for a known (AR t , A0 t , 
A 4>^ ) , using Equation (3.104). To obtain a diffraction point (V g , V p ) 
for a given receiver location (R^, 9^, <^), one can always assume the 
first diffraction point is at the source (V fi , V p ) * (V es » 
radiation direction (9^,$^ = ^) ^ or the positive ray (in +y direction) 

° r 55 2 ^) for a negative ray (in -y direction), and gradually add 

the increments (AR^,A0^,A^) until the final radiation direction 
(9 <!> ) is reached. The detailed description on this process was 

w W 

already given for the elliptic cylinder perturbation in Section III-C.2. 
One need not start out from the source everytime, but obtain the new 
diffraction point directly from Equation (3.104), provided that the new 
receiver location does not deviate greatly from the previous direction. 
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After the geodesic path is determined by this procedure, various 
other parameters associated with actual field calculations must be 
found. The Fock parameter 5 is given by 


5 « 



d £ 


(3.120) 


From Figure 3.10(b), the geodesic arc length is obtained as follows: 


i cosy * (r - S g ) sina . 


(3.121) 


or 


l siny = r $ - (r - S g ) cosa . 


(3.122) 


Therefore, 


r s cosy 

= cos^(Y-a) da * 


1 dS e 

~ sin(Y-a) d ^e * d ^e 


(3.123) 


(3.124) 
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The Fock parameter 5 is obtained by integrating along V r or V e , i.e.. 


a l^p 1/3 x 

5 - r s cost / UJ (-/) cos'll r-a') <lt ‘' 

V r 1 kfg 1/3 1 da 

' r s cos * ,j ^ (— ) cos^(r-a) dv: dv r 

v rs 


where 


. 1/2 

da_ 

dV r a 2 cos 2 V^ + b 2 sin 2 V p + z 2 cot 4 V e$ 


(3.125) 


(3.126) 


or 

V s x kpq l /3 1 dS e 

5 ■ / TJ ( 2 ) sln(-r-a) dT dV e * <’•«*> 

v es 

Note that p g * l/(<iCos2g + i^sin^B) and Kj, K 2 are two principal 
curvatures given in Equations (2.79) and (2.80), respectively. 

This completes the elliptic cone perturbation solution for the 
antenna mounted off the mid section of the ellipsoid. 
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4. Comparison with Exact Solution 

Geodesics for a source mounted on an ellipsoid can be analyzed 
precisely by computing the geodesic path defined by the surface 
parameters ( 9 q»$q) ar >d the geodesic tangent defined by the radial vector 
direction (9 t ,<j> t ) as shown in Figures 3.13(a) and (b), respectively. 

The geodesic path indicates the actual diffraction point location on an 
ellipsoid; whereas, the geodesic tangent indicates the radiation 
direction at the corresponding diffraction point. 

In order to show the validity of the perturbation solutions, a 
more elaborate numerical method for the geodesics employing calculus of 
variations, whose results are indicated as exact solutions here, is also 
studied and derived in Appendix A. Although this method provides 
accurate geodesics on the ellipsoid, it is far too complicated and 
inefficient to use for practical radiation applications. However, the 
exact solution is most appropriate for coupling problems where the exact 
path is desired between two known points on the surface. 

To show the validity of the elliptic cylinder perturbation 
solution, the source is placed at 0 S * 90° an d the geodesic paths and 
geodesic tangents associated with this source location are calculated as 
shown in Figures 3.14(a) and (b), respectively. In each figure, the 
elliptic cylinder perturbation solutions are compared with the exact 
solutions. Note that the 4X x 6X x 40X ellipsoid was chosen for more 
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(a) Geodesic path defined by the surface parameters ( 9 q* ^q) ■ 



(b) Geodesic tangent defined by the radial vector direction 
( » ♦t ) • 


Figure 3.13. Definition of geodesic path and geodesic tangent. 
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PERTURBATION 

EXACT 



(b) Geodesic tangents defined by the radial vector direction 

(v*t>. 


Figure 3.14. (Continued). 
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realistic aircraft and missile shapes. It is also noted that y is the 

A * 

angle between the geodesic tangent t and the principal direction t^ at 
the source location. Looking through those figures, one can see that 
the geodesic paths and the geodesic tangents of both methods coincide 
with each other within the significant energy region close to the 
source. This coincidence in the significant region can be checked more 
precisely by calculating the Fock parameter ( 5 ) along each geodesic 
path. Actually, the Fock functions associated with the solutions drop 
8.6 dB, 16 dB, and 23 dB as the Fock parameter (C) reaches 1, 2, and 3, 
respectively, in the deep shadow region. This clearly shows the 
significant portion of the surface as discussed in Section III-A. 
Although one can see small discrepancies for the rays toward the tips of 
the ellipsoid, i.e., lines for y = 80° and y = 280° in these figures, 
they only happen when the caustic effects come into play. In the 
caustic region where virtually an infinite set of rays have significant 
effects on the radiation pattern, the basic GTD theory fails. The study 
of this caustic effect is beyond the scope of this study. If one 
neglects the caustic regions, the geodesic paths and geodesic tangents 
using the perturbation models coincide with the exact solution very well 
in the significant region. 

The elliptic cone perturbation solutions can be examined by 
placing the source at 05*30° as shown in Figures 3.15(a) and (b), where 
the geodesic paths and geodesic tangents of the elliptic cone 
perturbation solutions are compared with those of the exact solution 
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(b) 


Geodesic tangents defined by the radial vector direction 
( e t t<J>t)* 


Figure 3.15. (Continued). 
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for the same 4X x 6X x 40X ellipsoid model. One should note the good 
agreement between the two results. The significant energy region of the 
geodesic path is also shown using the calculated Fock parameters (?) in 
the figures. For more detailed analysis with various source locations, 
one is referred to references [13,14], 

To make sure that the solution switches correctly between the 
elliptic cylinder and the elliptic cone perturbation method, the source 
is placed at 9 S =85 0 and the geodesic paths are calculated as shown in 
Figure 3.16. Comparing Figures 3.16 with 3.14(a), one can see that the 
geodesic paths of the elliptic cone perturbation solution are very close 
to those of the elliptic cylinder perturbation solution for that source 
location. 

These comparisons illustrate that the geodesic paths can be solved 
by using either numerical technique; however, the perturbation is much 
more efficient. In addition, one can easily relate the radiation 
direction with the desired geodesic path using the perturbation method. 
On the other hand, one is not sure which geodesic is necessary to 
achieve the desired radiation direction using the exact solution. 
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Figure 3.16. Geodesic paths defined by the surface parameters (e Q ,<f>Q) 
for a source mounted at 9s=85° on a 4X x 6X x 40X 
ellipsoid. 
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D. FLAT PLATE ATTACHMENT TO AN ELLIPSOID 


The flat plates used to simulate the aircraft wings and stabilizers 
are described by defining the location of the plate corners. The 
corners are numbered in a counter clockwise sense with the first and 
last corners being on the ellipsoid. The attachment of the flat plate 
to the ellipsoid is essential in that the intersection represents the 
wing-root section of the aircraft. The plates can be attached to the 
ellipsoid as illustrated in Figure 3.17(a) and (b). However, when the 
plates are attached on the lower half of the ellipsoid, the y component 
of the first and last corners are set equal to the y dimension of the 
ellipsoid center line as shown in Figure 3.17(c) and (d). This 
modification is made to eliminate the interior wedge problem which leads 
to multiple reflections and diffractions. Since the antenna is 
restricted to be on the top or near the top of the ellipsoid, this 
modification is reasonable. In addition, the field contribution from 
this edge is small compared with the other edges of the plate. This is 
due to the great attenuation of the surface wave propagating along the 
ellipsoid surface in reaching the edge. Note that the plates are not 
restricted to lie horizontally. The flat plate attachment above and 
below the center line of the ellipsoid will be discussed separately in 
this section. 

Now, the first approach to this problem is to find the intersection 
point between a line (i.e., one edge of the plate) and the ellipsoid. 
Then one can follow the same concept to find the curved junction edge 
between a flat plate and the ellipsoid. 
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Figure 3.17. Fuselage and wing geometries for aircraft model looking 
from the front. The antenna is assumed to be on the top 
portion of the models. 
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Using the geometry as showo in Figure 3.18(a), the ellipsoid 
surface is defined by 

R( v e , v r ) * a cosV e cosV p x + b cosV g sinV^ + c sinV fi z (3.128) 

and the unit edge vector is given by 

•\ ■ 7fy7 ■ 5 e lx + * e ly + 1 e iz <3 ‘ 129) 

I P 2" K 1 1 

where P^.y^Zj) and P 2 (x 2 ,y 2 ,z 2 ) are the position vectors of two 
corners with respect to the origin of the coordinate system. From 
Figure 3.18(a), the position vector of the first intersection point P HI 

can be easily defined by 

P„I • ?, - .ft . (3.130) 

Since pHI^ x h ,y h ,z h^ is on the surface of an el-|i P soid » one obtains the 
following equations: 

a cos V g cosV p = x x - tje lx 
b cosV e sinV p * y x - tje ly 

and 

c sinV e = Zj - tje lz . (3.131) 

From these equations, one finds 
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(b) Geometry of the curved junction edge. 
Figure 3.18. Flat plate attachment to an ellipsoid. 
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<T>> 



( x r*i e ix) 2 ( y i“ l i e iy )2 (z r £ i e iz^ 


(3.132) 


Defining 


A = b 2 c 2 e 2 x + a 2 c 2 e 2 y «■ aVe 2 z 


iy 


B = - b 2 c 2 Xl e lx - aVy^- a^z^ 


2 . 2 . 


.2u2. 


(3.133) 


C * b 2 c 2 x 2 + a 2 c 2 y 2 + a 2 b 2 z 2 - a 2 b 2 c 2 

and employing Equations (3.132) and (3.133), one obtains the distance 
between the points Pj and Pm such that 


*1 = 


-B + l B 2 -AC 
A 


(3.134) 


and 



The smaller one of H v t z ) will be used in Equation (3.130) to define 
the position vector ? HI (x h ,y h ,z h ). One can follow the same procedure to 
obtain the second attachment point P Hp . Then, as shown in Figure 
3.18(b), one can define the unit vector 


e Fl(x,y,z) 


P HF (x hf ,y hf ,z hf ) - P H I (x hi» y hi» z hi ) 

+ ♦ . . 
l P HF^ x hf » y hf * z hf ^ " p HI (x hi* y hi* z hi ) ' 


(3.136) 
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and the binormal unit vector which is given by 


t b (x,y,z) = e FI (x,y,z) x N p (x,y,z) 


(3.137) 


where the unit normal vector of the plate is given by 


e 1 (x 1 ,y 1 ,z 1 ) x e 2 (x 2 ,y 2 ,z 2 ) 

A 

N D (x,y ,z) * T I 

P |e 1 (x 1 »y 1 ,z 1 ) x e 2 (x 2 ,y 2 ,z 2 )| 


(3.138) 


According to the variation of the ellipsoid surface, one can divide the 

I 

line into N-l unequal length segments (t^) with i * 1,2...,N-1. 


The position vectors P^ along P HI P H p are given by the recursive equation 
such that 


P .j ~ P.j _ ^ + 4^ 6pi» i a 2,...,N-1 


(3.139) 


By using the position vectors Pj just found, one can attain N position 
vectors along the curved junction edge,1.e.. 


Pq£^ * P^ + 4^ t b for 1 3 2,3,....,N-1 


(3.140) 


with 


p CEi 


HI 


and 


PCE N 


HF 


(3.141) 
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Note that can be found by using the same idea employed to compute * r 
Therefore, the curved junction edge is found by connecting all N 
position vectors P(^(x,y,z). 

Next, consider the plate attachment below the center line of the 
ellipsoid. The first intersecting point Pm( x h* y h ,z h) can be defined 
by 



with 


x h “ X 1 ‘ l l e lx 


y h = y l ' Z I e ly 


and 

z h = Z 1 " *1 e lz 


(3.143) 


Considering the attachment to the elliptic cylinder as shown in 
Figure 3.19, one finds 

( y r*i e iy) 2 ( z r*i e iz) 2 

—i— + ;; » 1 . (3..144) 
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Defining 


A 

B 

C 


■ ♦ »Mz 

o ? 

= -c v,e,.. - b z,e 


y l e ly ' u "l c lz 
> 

1 • ~' z l 


c 2 y 2 + b 2 z 2 - b 2 c 2 


(3.145) 


and employing Equations (3,144) and (3.145), one obtains the distance 

-► 

between the points P^ and P^j such that 


-B + J B Z -AC 
4 1 A 


(3.146) 


and 



(3.147) 


The smaller one of U lt t 2 ) is used in E<l uation (3.142) to define the 

position vector P^i^h^h* 2 }))* ® ne can ^°^ ow same P rocedure 
obtain the second attachment point P^p* The procedure to obtain the 
curved junction edge is the same as the attachment above the center line 
of the ellipsoid. 

This completes the derivation of the curved edge junction between 
a flat plate and the ellipsoid. Note that the plates simulate the 
wings, stabilizers, etc. 
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CHAPTER IV 


RADIATION FIELD CALCULATION ALGORITHMS 


A. INTRODUCTION 

This chapter presents details of the techniques to determine each 
of the individual Uniform Geometrical Theory of Diffraction (UTD) terms 
which contributes to the total field. It is noted that, since UTD is a 
high frequency method, the lower frequency limit of this solution is 
dictated by the electrical dimensions of the ellipsoid and plate 
structures, i.e., the semi -mi nor axes of the ellipsoid are required to 
be at least a wavelength and each plate should have edges at least a 
wavelength long. In addition, each antenna element should be at least a 
wavelength away from all edges. With the antenna mounted on the 
fuselage it is obvious that the rest of the structure is very close to 
the soruce. That is to say, the antenna is in the near field of the 
structure. However, it is assumed that any point on the scattering body 
is in the far field of the source. 
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As depicted in Figure 4.1, the total radiated field is the 
superposition of the following UTD field components: 1) direct field 

from the source; 2) curved surface diffracted fields from the composite 
ellipsoid (referred to as the source field here); 3) reflected fields 
from the finite flat plates (i.e., wings or stabilizers); 4) diffracted 
fields from the edges of the plates, including diffraction from the 
curved junction edges formed by the intersection of the plates with the 
composite ellipsoid; 5) and vertex diffraction from each of the plate 
corners. In addition to the above direct scattering from a structure, 
the rays reflected or diffracted from one structure tend to interact 
with the other structures causing various higher order UTD terms such as 
double reflected, reflected-diffracted, diffracted-reflected, and 
double diffracted, etc. Thus, the total radiated field of an antenna in 
a complex environment can be expressed as a summation of individual UTD 
terms as follows: 

E "total ' h * E r + E d * h * E rr + E rd + E dr * E dd • 

Superimposing all of these scattered field terms is not unduly time 
consuming in that only a few contribute significantly in a given 
radiation direction. 

The first step in obtaining the solution of each UTD term is to 
determine the ray path using the laws of reflection and/or diffraction. 
After the ray path is identified, one must then examine the total ray 
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path to see whether or not It intersects an obstacle. If the ray path 
is not interrupted, the field value is computed and superimposed with 
other terms. On the other hand, if the ray path is interrupted, the 
field is not computed. A complete computer simulation program has been 
written which takes advantage of this feature and was used to provide 
the calculated results presented in the next chapter. 

Each UTD mechanism, which is expressed in terms of the x-y-z 
components of the electric field, will be discussed in detail in the 
following sections. 

B. SOURCE FIELD 

This section presents details of the techniques employed in the 
computer routines to obtain the direct and curved surface diffracted 
fields for a near field or a far field receiver. Both fields will be 
referred to as the source fields in the rest of this dissertation. The 
surface is assumed to be perfectly conducting and the surrounding medium 
is free space. An exp (jut) time dependence will be assumed and 
suppressed in the following formulations. 

The solution of the radiated field for an infinitesimal electric 
or magnetic dipole (i.e., dP e (Q’) or dP m (Q'), respectively) mounted on a 
perfectly conducting convex surface was presented in Sections II-B.(1,2) 
and the pattern factors P e and P m for an extended aperture (i.e., slot) 
and linear antenna (i.e., monopole), respectively, were given in Section 
II-B.3. Thus, the source field solutions can be obtained by replacing 
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dP e (Q') and dP m (Q') with the pattern factors P g and P m , respectively. 
The source field solutions are then summarized as follows: 


a) when the receiver is in the lit region. 


(i) for aperture type sources. 


,-jks 


E* = (L"n + L^b) + . 
m m m | jkp p 


in the near field, and 


‘m 


s r in the far field 


(4.2) 


where 


-jk 


L m = ^F [(P m *b’)(H £ + TqP cose 1 ) + (P m *T ' )T q F cose 1 ] , (4.3) 



(4.4) 


(ii) for monopole type sources. 


K = (Ljn + Lgb) 


-jks 


in the near field, and 


jkP »P 

e s r in the far field 


(4.5) 


122 



(4.6) 


-JkZo - 


"e 4 it e 


P. sine 1 [H^+T^Fcose 1 ] 


•JkZ 0 . 1 

~a P* sine 1 ! F 

4ir e o 


b) When the receiver is in the shadow region. 


(i) for aperture type sources. 


E* = (D"n + D^b) 
m m m 


,-jks 


/s(p c +s) 


in the near field, and 


in the far field 


- 1/6 




P g(Q‘) 

W 


- 1/6 


D m = + 


-jkt 


P g(Q') 

■POTT 


(4.7) 


(4.8) 


(4.9) 


(4.10) 
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(11) for the monopole type sources. 


(Djn + Dgb) 


e -jks 

/ / ^ i In the near field, and 

i s i p c +s ) 

j k P es *P r 

e es r in the far field 


(4.11) 


- 1/6 


°e- 


-jkZ 0 - 

4ir • e 


P He“ jkt 


P g(QV) 

TOTT 


(4.12) 


- 1/6 


»s- 




C g(Q') 


(4.13) 


It is noted that the caustic distance, p c , [30] is approximated by the 
geodesic arc length given in Equations (3.71), (3.121) and (3.122). It 
is a reasonable approximation in that the dimensions of the ellipsoid 
are assumed to be large in terms of the wavelength. The other 
parameters used in the above equations were defined in Section II-B. 

Through an extensive study [7] of geodesic paths on a surface of 
revolution, the dominant rays needed to be considered would not exceed 
four rays. These four rays are illustrated in Figure 4.2. However, the 
significant effect of the ellipsoid surface is associated with a region 
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1 ^"LOCATION 


DIFFRACTION 
POINTS 



Figure 4.2. The four dominant 6TD terms that radiate at 
(9=90°, ♦-145°). 

around the source. Since the fuselage of the usual aircraft is very 
long and slender, only two dominant rays are included in the 
computation of the source field in the deep shadow region. It- is also 
observed that rays 3 and 4 in the figure fall into the caustic region 
where the basic UTD theory fails. In order to avoid this caustic 
problem, the cone boundary shown in Figure 4.3 is used in determining 
whether one or two rays are used in the solution. Note that Sj2 
defined automatically by determining the caustic angle In the elevation 
pattern (8 C ) and adding a few additional degrees to that value, i.e., 
g 12 * e c + 48 where 2° < 4S < 10°. The caustic angle S c can be 
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LIT REGION 


TANGENT PLANE 



Figure 4.3. Cone boundary used to define terms to be included in 
the shadow region. 

obtained by enforcing the determinant of the first matrix in Equation 
(3.104) to zero with the following assumptions: 

r s - S e 0£ » y * S * ir/ 2 • a • 0 and ^ » 0 or t • 

Note that the above parameters were defined in Section III-C.3. One 
would expect to observe slight discontinuities somewhere, because 
various numbers of rays are included in different regions. 

This completes the source field calculation and this source field 
will also be used as an incident field at the scattering body of the 
aircraft structure in the following sections. 
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C. REFLECTED FIELD 


The reflected field from each of the finite flat plates is 
determined using image theory [24], As shown in Figure 4.4, the 
receiver image position with respect to plate is given by 





(4.14) 


The source field solution is, then, used to compute the electric field 

♦I . 

(E S ,E S ,E S ) at the image position P p J . This field would exist at that 
x y z 

point if the plate were not present. The second step is to ascertain 
if the ray path from the effective source to the image position 
intersects the plate. This ray path shadowing algorithm will be 
discussed in detail in Section IV-F. If the ray path is not 
interrupted by the plate, a reflection by that plate does not occur 
for that receiver position. If the ray does pass through the plate, 
then the reflected field (E r = E^x + E^y + E^z) is given by 


E x< f r> 


T xx T xy T xz 


E x<% j > 

e;oy) 

X 

T yx T yy T yz 


eJ(?J J > 

E ;< ? r> 


T zx T zy T zz 

; 

e|(p’ j ) 

1 

1 


— 


— — 


(4.15) 
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Figure 4.4. Geometry used to determine reflected field from 
plates (wings, etc.). 



where [T] represents the reflected field polarization transformation 
matrix which satisfies the boundary conditions [32] on the given plate. 

A A 

This matrix is determined using the normal n and tangent t unit vectors 
of the electric field to the given plate as follows: 


E r * (n • E s )n - (t • E )t on p] a t e (4.16) 
(t • E s )t = E s - (n • E s )n (4.17) 
E r = 2(n • E s )n* - I s (4.18) 


or 


A A A A A A A 

T xx x + T^ x y + T zx z * 2(n • x) n - x 

A A A A A A A 

T X y x + T yy y + T zy z ■ 2(n • y) n - y 



(4.19) 

(4.20) 

(4.21) 


Since the T-matrix is independent of the receiver location, it is stored 
for each plate in order to optimize the computational efficiency of the 
numerical solution. The total reflected field for a given receiver 
position Is, then, the superposition of the reflected fields from the 
individual plates. 
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D. DIFFRACTED FIELD 


The diffracted field is obtained by using the UTD formulation 
presented in Section II-C. The shadowing of the incident and reflected 
rays by the flat plates is compensated by the addition of the diffracted 
fields which cause the total field to be smooth and continuous at the 
incident and reflection shadow boundaries. The most difficult part in 
computing the edge diffracted field is to determine the ray path. For 
any particular source, plate edge, and receiver location, the diffracted 
ray path is unique [25], The key in finding this unique path is to 
determine the actual diffraction point along the given edge. 

Typical geometries for the diffraction point in the illuminated and 
shadow region are depicted in Figures 4.5(a) and (b), respectively, and 
will be discussed separately in the following sections. 

1. Diffraction Point in the Lit Region 

As shown in Figure 4.5(a), the diffraction point in the lit region 
must be located at a position along the edge such that the angle (3 
between the incident ray and the edge equals the angle 8* between the 
diffracted ray and the edge. With i and d defined as unit vectors in 

A 

the direction of the incident and diffracted rays, respectively, and e 
defined as a unit vector lying along the edge, the above noted basic law 
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(b) Diffraction point in the shadow region. 

Figure 4.5. Illustrations of the diffracted field by the straight 
edge. 
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of edge diffraction [27] (i.e., 6 - 3’) can be stated equivalently by 
the requirement that 

i • e = d • e . ' (4.22) 


The edge unit vector e in Equation (4.22) can readily be computed from 
the specified corner coordinates of the plate. From the procedure 
depicted in Figure 4.6, the diffraction point is given by 



4 









in the near field case. 


(4.23) 


or. 


4 4-1 



4 


+ IP, 


4-1 a 

- P | cose e i 


in the far field case. 


(4.24) 


where 

4-1 4 * 4 - 4 - * * 

p s » p Ci ♦ {(p s - p c .) • «,} 

P r - K + {(P r - P c ) • e i) e i 

r i i 

and 


(4.25) 

(4.26) 


cote * 




(4.27) 
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2. Diffraction Point in the Shadow Region 


If the edge of the plate is not located in the lit region of the 
source as shown In Figure 4.5(b); the diffraction point cannot be 
obtained directly in that one has to resort to a numerical approach. An 
iterative numerical search procedure is employed here to determine the 

A A A A 

diffraction point which satisfies the diffraction law (i.e., i*e = d*e). 
As an initial step in the routine prior to initiating the actual search 

A A 

procedure, values of i »e are computed and stored for a sequence of N 
sample points along the plate edge. The coordinates (x,y,z) of the 
sample points are selected and defined in the following manner. For an 
edge of length £, the N sample points effectively subdivide the edge 
into (N-l) segments of length A£, where A£=£/(N-1). The distance t 
between the one reference end of the edge and the n^ sample point along 
the edge is then given by 

t(n) = (n-l) A£ ; n«l,2,3 ,N . (4.28) 

These distances define the coordinates (x n ,y n ,z n ) ^ or eac ^ t ^ ie N 

A A 

sample points. It is noted that the i«e curve for the straight edges 
which attach to the fuselage varies rapidly near the fuselage. So for 
these edges, the sample points are taken in a nonlinear way with more 
samples near the fuselage, i.e., the distance t is given by 

n-l 

t(n) * J i As j n a l , 2,3 , • • • . ,N (4.29) 

i*0 
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where 


t 

AS = . ( 4 . 30 ) 

l i 

i*0 

Except for this variation, the routine follows the detailed approach 
which was given previously. 

After the i»e values are computed and stored at the N sample points 
along each plate edge, the routine initiates a search to determine which 
one (if any) of the (N-l) segments contains the diffraction point. The 
test to determine whether the diffraction point is within the limits of 

A A 

the physical edge or not can be done by computing d»e at two ends of 
each physical edge for a specified receiver point and then comparing 
them with stored data i *e at the two ends. A diffraction point lies 
along the physical edge if d # e>i*e at one end of the physical edge and 
d*e<?*e at the other end. Otherwise, the diffraction point does not 
exist on the physical edge. If the diffraction point falls within the 
limits of the physical edge, the routine selects the midpoint of the 
edge as the next test point. Note that a selected test point in the 
search routine always coincides with one of the N sample points in the 
list of previously computed i «e values. Noting the sign of the 
inequality between i »e and d»e at the midpoint and two endpoints of the 
edge, the routine selects the half of the total edge which contains the 
diffraction point as the second search interval. The routine continues 
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to halve the search interval in this manner until, after several cycles 
in the iteration, the length of the search interval is reduced to the 
segment length At between successive sample points. This situation is 
depicted in Figure 4.7, which shows the diffraction point located within 
the segment bounded by distances t(m) and t(m+l). At this stage of the 
search, the routine uses linear interpolation to compute the distance 
h d to the true diffraction point. The value of h d is then used to 
compute the diffraction point coordinates P d (x d , y d , z d ), as given by 


x d ■ X 1 * h d • e x 
U * h * h d * e y 


and 


z d = 


Z 1 + 


(4.31) 

(4.32) 

(4.33) 


3. Curved Edge Diffraction Point 

This section presents a discussion of the routine for computing the 
diffraction point on a curved edge of a finite flat plate. In general, 
a curved edge solution is needed in the analytical aircraft model to 
treat the diffraction which occurs at the junction of the finite flat 
plates and the ellipsoid, and the chopped fuselage curved edge. 
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The basic approach used to determine the diffraction point along 
the curved junction edge is essentially the same as that for the 
straight edge diffraction point i/i the shadow region which was discussed 
in the previous section. The only difference is that the edge unit 
vector e is a function of position along the curved edge and the 
determination of the incident unit vector i is more complicated in that 
the curved junction edge lies on the ellipsoid surface as shown in 
Figure 4.8(a). The position vectors of the N sample points along the 
curved junction edge are easily computed using the attachment algorithm 
which was discussed in Section III-D. Accordingly, the i th edge unit 
is obtained as follows 





(4.34) 


The unit vector i is determined in the following manner. Since the 
curved junction edge lies on the ellipsoid surface, one must use the 
elliptic perturbed cylinder or elliptic cone method which was discussed 

a 

in Section III-C to find the incident unit vectors. The unit vector i 
for each of N sample points along the curved junction edge is expressed 

by 


(4.35) 

i * t : cose + t e sine 


/ 
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(b) Incident unit vector on the curved junction edge. 



(c) Incident unit vector on a developed elliptic cone. 
Figure 4.8. (Continued). 
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where 0 denotes the angle between t^ and the unit vector i. The 
location of the source and the i^ sample point on the curved junction 
edge are defined as P s (V es ,V ps ) and respectively. 

For the antennas mounted on the off-mid section of an ellipsoid, 
the perturbed cone method is used in the calculation of incident unit 
vectors as shown in Figures 4.8(b) and (c). The 0 angle is, then, 
given by 

B - y - a , ( 4 *- 

where a was given in Equation (3.86). The y is obtained from the 
following geodesic equation: 


r cos(y - a) * r cosy . 


(4.37) 


where 


= (a* cos 2 V + b 2 sin 2 V + z 2 cot 4 V ) 1 ^ , 


S ' s 


rs s rs s , *''“ es- 


(4.38) 


r e = (a 2 cos 2 V ri + b 2 sin 2 V p1 + z 2 cot 4 V es ) W2 - , 


(4.39) 


and 


/ e1 [c 2 cos 2 Vg + (a 2 cos 2 V pi + b 2 sin 2 V r1 )sin 2 V e ] 1/2 » dV e , 


es 


(4.40) 


with 
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and 


a s * a cosV es 
b $ * b cosV es 


z $ * a sinV e$ 


(4.41) 


Thus, y is given by 


-1 / s - r e cos ° 
tan l r P sina J 


(4.42) 


A 

Combining Equations (4.35), (4.36) and (4.42), the unit vector i at 

each of N sample points along the curved junction edge is then computed. 

For the antennas mounted on the mid-section of an ellipsoid, the 
elliptic cylinder perturbation method is used in the calculation of the 
incident unit vector i along the curved edge. Note that the angles 6 
and y are the same for the elliptic cylinder case. The geodesic 
equation is, then, given by 


S e 


B = tan 


(4.43) 


where 


* / el {c 2 cos 2 v‘ + (a 2 cos 2 V r + b 2 sin 2 V p )sin 2 V e } 


o 9 o •, 1/2 * 


(4.44) 


142 



and 


S r * / 

r V 


H 


cosV 


es 


rs 



(4.45) 


Combining Equations (4.35) and (4.43), the incident unit vector i at 
each of the N sample points along the curved junction edge is, then, 
computed. 

After the edge unit vector and the incident unit vector at each of 
the N sample points on the curved edge are computed and stored, the 
actual diffraction point on the curved edge is determined using the 
diffraction point search algorithm for the shadow region. 


4. Diffracted Field Calculation 


Assuming that the ray path is identified, the diffracted field for 
the m th edge of the j th plate is given by 


E 1< ? r> 

E 1<V 


D s (L,* , ,*,B 0 ,n) 0 

0 D h (L,*',<M 0 ,n) 



The parameters needed for the diffracted field calculation are given as 
follows: 


s d (s^ + t) 
L * s d + s 1 + l 


sinV 


(4.47) 
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B 0 = sin -1 ( 


(4.48) 



(4.49) 


(4.50) 


(4.51) 


and 

s d in the near field, and 

-Pp'Pd in the ^ ar (4.52) 

Note that i s (? d ) is the field incident on the edge at the diffraction 
point and is obtained using the source field solution which was 
discussed in Section IV-B. The incident field can then be expressed as 

E S ,(P d ) - E n (n-e*) + E b (b-{') (4.53) 

and 
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The unit vectors of the edge fixecj coordinate system for the incident 
ray are given by 


♦ ' = 


n. cos 4>* 

J 


- b„ si n <(>' 


(4.55) 


and 


B‘ = *' x 1 . (4.56) 

Now, the diffracted field in terms of the rectangular coordinate system 
for the m^ edge of the plate is given by 




♦ ♦*1 


(4.57) 


with 

A ^ ^ 

* * n, cos* - b_ ^ sin* (4.58) 

J 

and 

8 - * x d . (4.59) 

These unit vectors (*',6',*,S) will be used in the following sections 
without defining them again. 


5. Corner Diffracted Field Calculation 


The corner diffracted field is usually very small except near the 
edge diffraction shadow boundaries. Therefore, the physical edge length 
is extended by 5 wavelengths for the corner diffracted field calculation 
in that the current corner diffraction coefficients are based on the 
concept of equivalent edge current [31,333, which would exist on this 
edge when extended to infinity, as discussed in Section II-C.2. 

There are two corner diffraction terms associated with each finite 
plate edge. However, for the edges attached to the fuselage, there is 
only one corner diffraction term involved in the calculation. The 
corner diffracted field and its related parameters are defined in 
Section II-C.2. Some parameters related to the geometry given in Figure 
4.9 are defined as follows: 

*c = cos ’ 1 ( *c * e m,j ) * ( 4 * 6 °) 

6 oc = cos_1 (s * * e m,j } < 4 * 61 > 


and 

s c s 

* Sc + ^ * (4,62) 

The total corner diffracted field is obtained by superimposing the 
corner diffracted fields from each of the corners. 
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Figure 4.9. Illustration of the corner diffracted field. 
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E. FIELD CALCULATION FOR HIGHER ORDER TERMS 


In addition to the above direct scattering from a structure, the 
rays reflected or diffracted from one structure tend to interact with 
the other structures causing various higher order UTD terms (i.e., 
double reflected, reflected-diffracted, diffracted-reflected, and 
double diffracted). 

1. Double Reflected Field 

The double reflected field E rr can be obtained following a 
procedure similar to the single reflected field in Section IV-C. As 
shown in Figure 4.10, the images of the receiver location with respect 
to the j th plate and j th -k th plates are given, respectively, by 

P* j = K “ 2 "j * (P r ' ^c • (4.63) 

1 * J 

and 

PpJ.K = p^J _ 2r , k . (p^J . P c ) . (4.64) 

1 ,k 

The source field solution is, then, used to compute the electric field 
E s at the image position P* 3,k . This field would exist at that point if 
the j th and k th plates were not present. Then the double reflected 
field is, then, given by 

[E V V r )l ■ [T]j [Tl k [i s ($' k )1 • (4.65) 
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Figure 4.10. Double reflection geometry 



Note that [T]j and CTD^ satisfy the boundary conditions on the j th and 
k^ plates, respectively, and the matrix [T] was defined in Section 
IV-C. 


2. Reflected - Diffracted Field 

A typical geometry for the reflected diffracted field is 
illustrated in Figure 4.11. To begin the solution, one needs to define 
the images of the second plate corners and the receiver location with 
respect to the first plate. From the geometry of Figure 4.11, the image 
point of the receiver with respect to the k th plate is given by 


Pj k = P r - 2n\ • (P r - P Ci i( ) . (4.66) 

The images of the m th corner of the m th edge of the j th plate with 
respect to the k^* 1 plate are given by 






and 


-Ik 


e m,j 


P 


IP 


Ik 

c m+i , j 
Tk 

C m+i , j 



(4.67) 


(4.68) 
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respectively. The value i*e^ k i is then computed and stored at the N 
sample points along each edge of the image plate. The iterative 
numerical technique which was dis-cussed in Section IV-D is used to find 
the diffraction point P d k on the edge of the image plate. From the 
diffraction point on the image plate, the actual diffraction point at 
the m th edge of the plate is obtained as follows: 

+Il * +Ib + \ 

p d * p d + 2 "k • ( p d - p c lik ) • < 4 - 69 > 

After the ray path is identified, the source field E s is computed at the 
+j. -r^ 

diffraction point P d on the image plate. Next, the reflected field E 

at the actual diffraction point P d is obtained using the single 

reflection solution which was discussed in Section IV-C, and is given 

by 


[i rk (P d )l - [T] k 




(4.70) 


Then the reflected-diffracted field is given by 


” r,d • + ~ 
E, m ’ J (P p ) 


D s (L,*',+,B 0 ,n) 0 


r k * 
e , < p d> 

E ^k d m,j ( ^ r ) 

S m. 

0 (L » ♦' »$,B 0 ,n) 

1 

r k - 
E i< p d> 

— ^ — 


— — 


— — 



(4.71) 
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The parameters needed for the diffracted field calculation are given as 


follows: 


s rd (s r + s 1 + l) 
L “ s rd * s r + S 1 + i 


sin^B, 


(4.72) 


<)>' * tan 


-1 


A r * 
-s • 


-s 


r . b 


m,3 


(4.73) 


4 = tan 


-1 


:rd : 
s • n-i 

s rd - b " 


(4.74) 


B 0 - sin " 1 ( Jl- |s Pd • e m ,j | 2 ) • 


(4.75) 


A = 


s r ♦ s 1 


s rd (s rd + s r ♦ s 1 ) 


(4.76) 


and 


s * 


s rc * in the near field, and 


-P . p, in the far field 
r d 


(4.77) 
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The reflected-corner diffracted field can be computed following the 
same procedure as that of the reflected-diffracted field except using 
the corner diffraction coefficient rather than the edge diffraction 
coefficient, and is not repeated here. 

3. Diffracted-Reflected Field 

The geometry for the di ffracted-reflected calculation is 
illustrated in Figure 4.12. To begin the solution, one needs to define 
the image point of the receiver location with respect to the j th plate 
which is given by 




(4.78) 


♦I,- 

If one assumes the image point P p J as the receiver location, the 
diffraction point along the m^ edge of the plate can be found 
through the single diffraction point search technique which was 

-Ii 

discussed in Section IV-D. The diffracted field at the image point P p J 
is, then, given by 



Dg(L»$' »+»8q,o) 0 



0 D h (L,t',<M 0 ,n) 








Ae 


-jks 


(4.79) 
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Figure 4.12. Diffracted-reflected geometry. 
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Some interesting parameters for this diffracted field calculation are 
given as follows: 



( 4 . 80 ) 


( 4 . 81 ) 


( 4 . 82 ) 


( 4 . 83 ) 


( 4 . 84 ) 


( 4 . 85 ) 
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Finally, the diffracted-reflected field is obtained by multiplying the 
reflected field polarization transformation matrix to the above 
diffracted field, and is given by. 


[ | d m.k r j(P r )] . [ T ] 




(4.86) 


For the corner diffracted-reflected field calculation, one needs to 
use the corner diffraction coefficient rather than the wedge diffraction 
coefficient. Beside that, the corner diffracted-reflected field 
can be computed following the same procedure as that for the 
diffracted-reflected field, and is not repeated here. 

4. Double Diffracted Field 

The double diffracted field contribution is usually small compared 
with the other UTD terms. However, when one edge lies in the shadow 
boundary of the other edge or both diffracted edges are close to the 
source, a double diffracted field may be significant. A typical 
geometry for the double diffracted field calculation is illustrated in 
Figures 4.13(a) and (b). The key to solving for the double diffracted 
field is to determine the two diffraction points along the diffracting 
edges for specific source and receiver locations. 
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EFFECTIVE 



(a) Double diffracted field when the two diffraction points 
lie on the different plates. 

Figure 4.13. Illustration of the double diffracted field. 
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(b) Double diffracted field when both diffraction points 
lie on the same plate. 


Figure 4.13. (Continued). 
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When the two diffraction points lie on different plates as shown in 


Figure 4.13(a), one needs to first compute and store the values 

e . . s * at the N sample points, along the second edge e The N 

m,j 1 

sample points along the second edge are assumed as the receiver 

A 

locations and the first edge e n ^ as the diffracting edge. Then the 

value e • • s'** at each sample point on the second edge can be computed 
m,j 

using the usual single diffracted field algorithm and stored. In this 

-► 

procedure, the corresponding diffraction point P d ^ on The first 

diffraction edge for each sample point on the second edge is also 

determined using the diffraction point search technique which was 

discussed in Section IV-D. Next the diffraction point P d2 on 

the second edge is found using the diffraction law, i.e., 

e • s** 1 = e • • s* 2 . After the ray path is identified, the double 
m,j 

diffracted field is computed in the following manner. The incident 
field E dm *^ at the second diffraction point is given by 


E^(P d2 ) 

s _ 

l M L i**i’V e oi» n i ) 0 



E^tPdj) 

1 

1 

0 

D<j (Tj» ^ » n ^ ) 


E*(? dl ) 


a 

• A.e 

1 (4.87) 

where E s is the source field, which was discussed in Section IV-B at the 
first diffraction point P dr The parameters needed for this diffracted 
field calculation are given as follows: 
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di ^ 

s V + l) 


■ 1 '71 


sin 2 6, 


s *+ s' + i 


ol 


(4.88) 


6 ol " sin 


( ji- I* 1 - VjI 2 ) . 


(4.89) 


1 _i 

= tan 


■s • 


•s 1 • b 




(4.90) 


<t> * tan 


-1 


■»di » 

s n< 


»d i » 

s *• b 




(4.91) 


and 


A = 


s dl (s dl + s 1 ) 


(4.92) 


Finally, the double diffracted field at the receiver location is given 
by 


E Vj d n,k ( ? r , 

S M 

1 

Ds(L 2 ,* 2 **2 f ®o2 ,n 2 ^ ® 


E^ mJ (P d2 ) 

E d m,J d n.k ( ^ ) 


0 1 

°s ^ 1 "2 * ^2 * ^2 * ®o2 ,n 2^ 


E^ mJ (P d2 ) 


Age 


- jks d 2 


(4.93) 
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The parameters needed for this diffracted field calculation are given as 
follows: 



(4.94) 



(4.95) 


(4.96) 


(4.97) 



and 

S d 2 



in the near field, 
in the far field 


and 


(4.98) 


(4.99) 
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If both diffraction points lie on the same plate as shown in Figure 
4.13(b), the ray path should be determined in a different way. The 
first step is to project the receiver onto the edge which contains the 
second diffraction point which is given by 


>1 




+ 




“m 




"m.j 


(4.100) 


Next, project the receiver onto the plane which contains the plate 

.di 

itself and the second incident ray path s . Then, the projection of 
the receiver is given by 


•Ml 

P r 


-► I 

P r 


f I A 

l p r - P r l », 


m,j 


(4.101) 


-Ml 

Assuming that the receiver is located at P p , one can find the first 
diffraction point P^ using the single diffraction point search 

technique which was discussed in Section IV-D. After the first 

+ ...•*■ 
diffraction point P^ is determined, the second diffraction point P^ 

is given by 



(4.102) 


where 


- cos’ 1 (b 


m,o 


s'" 1 ) 


(4.103) 
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and 





(4.104) 


After both diffraction points are determined, the double diffracted 
field can be computed in the same way as done for the diffraction points 
on the different plates. 

It is noted that the parallel component of the field incident on 
the second edge is always zero for this case, and if the diffracting 
plate is a half plane (i.e., n=2) there are two perpendicular fields 
incident on the second edge as shown in Figure 4.14. 

The diffracted-corner diffracted, corner diffracted-diffracted and 
double corner diffracted terms can be computed using the same procedure 
as done for the double diffracted field, however, these terms are 
usually very small. 

F. SHADOWING ALGORITHM 

After the ray path is determined using the previous field 
calcualtion algorithms, one must then examine the total ray path to see 
whether or not it intersects an obstacle. If the ray path is not 
interrupted, the field value is computed and superimposed with other 
terms. On the other hand, if the ray path is interrupted, the field is 
not computed. 
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Pes 



Figure 4.14. Two perpendicular components of the field incident to the 
second edge. 
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First of all, it is necessary to decide whether the incident ray is 
toward the plate. This can be done by determining the value of the 
parameter which is given by 

C 1 * t n j * (P r " P es^ * f n j * " p es^ * (4.105) 

i 9 J 

If the value of Cj is less than or equal to zero, such as the point Pj 
in Figure 4.15, the ray cannot be interrupted by the plate. Otherwise, 
a next step is needed. From the geometry given in the figure, the 
intersection point between the incident ray and plate of the plate 
is given by 


= p + 
es 


C 2 < f r - ? 


es' 


(4.106) 


where 





(4.107) 


If the value of C 2 is greater than one, then the ray cannot reach to the 
j ^* 1 plate such as the point P 2 as shown in the figure. Otherwise, one 
more step is needed to decide whether the ray passes through the plate 
which is defined by corners. For this test, the parameter C 3 is used 
and given by 
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- IT 


( 4 . 108 ) 


Nj 

I I J 

m=l 



» 


where 


0 m = tan 
m 


-i 


"j • t( p c m - P H> * ( p c mt . - P H>] 
m, j m+1 9 j 

_ < P C. P H> • ( p c_„ P H> 

Ml-J l»' r A»J 


( 4 . 109 ) 


and Nj is the total number of corners used to define the plate. If 
the value of C3 is less than zero, then the ray is not shadowed by the 
plate such as the point P3 in the figure. Otherwise, the ray is 
shadowed by the plate such as P4 in the figure. 

All ray paths for each UTD term are tested with this shadowing 
algorithm to decide whether or not the UTD term is included in the total 
field. This completes the total field computation. 
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CHAPTER V 


ANALYSIS OF AIRBORNE ANTENNA RADIATION PATTERNS 
A. INTRODUCTION 

Using the fuselage simulation techniques and the field calculation 
algorithms presented earlier, the Ellipsoid-Model Aircraft Code [15] has 
been developed to compute and superimpose the various UTD terms for an 
arbitrary near or far field pattern. 

To demonstrate the capability of this new analytic solution, 
various airborne antenna radiation patterns are analyzed using the code. 
Some radiation patterns that have been previously computed with old 
programs (Elliptic cylinder-model Aircraft Code [40,45] or Prolate 
spheroid-Model Aircraft Code [24]) are also computed and presented here 
in that this new solution provides an improved result in most cases. 

The validity of this analytic solution is also verified by the 
comparisons between the calculated and measured data for the wide 
variety of shapes such as commercial, private and fighter aircraft, 
missiles, and spacecraft. Most of these comparisons are in terms of 
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actual aircraft simulations that have been used in the past for 
verification purposes. The solutions are presented in terms of the 
radiation patterns of the three basic elements (i.e., axial slot, 
circumferential slot, and monopole) which can be used to simulate an 
arbitrary fuselage mounted antenna if the aperture distribution is 
known. 


T n nrHflr t A avami na ami* _j. . . . _ 

*" bulimic TUI IUUD uumcai paucm LUlb , a Cdrtesian 

coordinate system (x,y,z) originally defining the simulation model is 
now rotated into a new system (x p ,y p , 2 p) as shown in Figure 5.1. Note 
that the new cartesian coordinates are found by first rotating about the 
z-axis an angle <t» c and then about the y-axis an angle 0c. The pattern 
is, then, taken in the (x p ,y p ,z p ) coordinate system with 0 p fixed. The 
coordinate systems that apply for the principal plane pattern cuts 
(i.e., roll, elevation, and azimuth patterns) are illustrated in Figure 
5.2. 


To begin any simulation of an aircraft, one needs to start with a 
set of scale model drawings. For our purposes, Jane's book [41] "All 
the World's Aircraft" is most appropriate in that it provides a 
reasonable line drawing for a wide variety of aircraft. 


B. COMMERCIAL AIRCRAFT 


This section presents various radiation pattern results for 
different antennas, locations and commercial aircraft. Measured 
radiation patterns were provided by NASA (Hampton, Virginia) for a 
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y p 


(a) ROLL PLANE PATTERN COORDINATES (6 *0® ) 

c * c 



(b) ELEVATION PLANE PATTERN COORDINATES (^=90°, <£ c =90®) 


Figure 5.2. Transformed coordinate systems for the 
principal plane pattern cuts. 
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(c) AZIMUTH PLANE PATTERN COORDINATES (0 C *9O°, <£ c = 0°) 


Figure 5.2. (Continued.) 
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Boeing 737 and KC-135 aircraft. Since most commercial aircraft are 
similar to the 737 and KC-135, one can assume similar results would 
apply for other such structures.. However, if the aircraft of interest 
deviates to a high degree from these geometries, one must verify his 
results before he can conclude their accuracy. Once that new 
configuration is verified, the class of geometries that can be treated 
increases accordingly. 

1. Boeing 737 Aircraft 

The Boeing 737 aircraft as shown in Figure 5.3 is most interesting 
in that extensive experimental work [42,43,44] is available. 

One of the major problems in the design of a safe, reliable and 
accurate approach and landing system known as the Microwave Landing 
System (MLS) is the location of the antennas on the aircraft structure 
in order to achieve the desired radiation coverage. For the Beoing 737 
aircraft. Stations 220, 250, and 305 on the top of the fuselage and 
Stations 222 and 950 on the bottom of the aircraft, as shown in Figure 
5.4 were proposed as test locations for the antenna installation. In 
order to determine the best location for the MLS application, the 
radiation patterns for each test location was analyzed and measured. 

Figure 5.5 illustrates the computer model which is used to simulate 
the Boeing 737 aircraft for the top mounted antennas. A composite 
ellipsoid (77" x 74" x 830" x 308") is chosen to simulate the fuselage 
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ANTENNA LOCATION 

A - STATION 220 
B- STATION 250 
C- STATION 305 
D - STATION 222 
E - STATION 950 



Figure 5.4. Test locations for the antenna installation on the 
Boeing 737 aircraft. 
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surface of the actual aircraft as accurately as possible, especially 
near the antenna location. It is assumed here that the radome is 
perfectly transparent. This leaves a short, blunt-looking nose section 
which extends out from the front of the aircraft. It was found that the 
nose section normally has little effect on the resulting pattern. Thus, 
the cockpit nose section for simplicity is simulated by a finite flat 
plate. The wings and vertical stabilizers are approximated by finite 
flat plates which results in a simple model for the aircraft. The 
significant features associated with the vertical stabilizer are the 
leading edge and thickness of the two flat plates used to approximate 
those features of the actual vertical stabilizer. The thickness is 
important in that it tends to shadow the direct field from the antenna 
for aft radiation directions. For the sake of efficiency, no horizontal 
stabilizers are included in this model. 

The radiation patterns for a X/4 monopole mounted at Station 220 
above the cockpit on a Boeing 737 aircraft are, then, calculated using 
the model just described. The three principal plane results are shown 
in Figures 5.6 to 5.8 and found to be in very good agreement with 
measurements. The experimental work was performed by the technical 
staff at NASA (Hampton, Virginia) using a l/ll th scale model of a Boeing 
737 aircraft. But it is noted that the measured results have some 
asymmetry in the patterns. This could be attributed to misalignment of 
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CALCULATED 

MEASURED 



(SCALE: EACH D I V I S I ON* 10DB) 


Figure 5.6. Elevation plane pattern of a X/4 monopole mounted at 
Station 220 on top of a Boeing 737 aircraft. 




CALCULATED 

MEASURED 


TOP 



BOTTOM 

(SCALE: EACH 0 1 VIS ION* 1 OOB) 


Figure 5.7. Roll plane pattern of a X/4 monopole mounted at 
Station 220 on top of a Boeing 737 aircraft. 
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CALCULATED 

MEASURED 



NOSE 



TAIL 

(SCALE: EACH DIVISION* 1 ODB) 


Figure 5.8. Azimuth plane pattern of a X/4 monopole mounted at 
Station 220 on top of a Boeing 737 aircraft. 
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the monopole with respect to the surface normal or the movement of the 
model due to shifting weight during the measurement. This misalignment 
of the monopole (approximately 3° tilted to the nose section from the 
actual surface normal) was detected by the calculated elevation plane 
pattern and various azimuth plane patterns. To compensate this 
misalignment in the comparison with measured patterns, the conical 
pattern axes were rotated 3° toward the nose section. The various 
azimuth plane patterns for this antenna location are computed and shown 
in Figures 5.9(a) to (g). In each case, the calculated results compare 
very favorably with the measurements. 

To evaluate the radiation coverage performance of the MLS antenna 
at Stations 250 and 305, the elevation plane patterns for a X/4 monopole 
mounted at these locations are obtained. The elevation plane patterns 
are used in that it is the most critical MLS sector. It is noted that 
the antenna at Station 250 is mounted 4" off the fuselage centerline 
because of the structural problems that exist on the actual aircraft 
fuselage. Both calculated and experimental results for Stations 250 and 
305 are presented in Figures 5.10 and 5.11, respectively. The results 
reveal good agreement between the theoretical predictions and scale 
model measurements. In addition, the results indicate that Station 220 
appears to be the best choice since It gives the best forward coverage 
which satisfies the MLS performance requirement. However, due to 
structural problems, the antenna of interest had to be moved to 
Station 250. 
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NOSE 



TAIL 

(b) 0 p = 6O* 

Figure 5.9. Azimuthal conical patterns of a X/4 monopole mounted at 
Station 220 on top of a Boeing 737 aircraft. 
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NOSE 



TAIL 

(d) 0 p = 80° 


Figure 5.9. (Continued). 


184 


CALCULATED 

MEASURED 


RIGHT 

HING 
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( e) 0 p = 100° 
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TAIL 

(f) 0 p = IIO e 

Figure 5.9. (Continued). 
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CALCULATED 

MEASURED 


NOSE 



TAIL 

(SCALE: EACH DIVISION- 10DB) 


(g) 0 p = l2O° 


Figure 5.9. (Continued). 
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BOTTOM 

(SCALE: EACH 0 1 V IS ION= 1 ODB) 


Figure 5.10. Elevation plane pattern of a X/4 monopole mounted at 

Station 250 (off center) on top of a Boeing 737 aircraft. 
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Next, an antenna location on the bottom of an aircraft fuselage 
behind the radome was also proposed for the MLS application. Since the 
MLS is used for both for the approach and landing, this location 
appeared to be very attractive for the forward coverage requirement. 

For an antenna location at Station 222, the Boeing 737 aircraft was 
modeled by a composite ellipsoid (66" x 55" x 785" x 232"). It is noted 
that the fuselage dimensions associated with the bottom of the aircraft 
are different from the top, in that the aircraft has a definite shape 
change from top to bottom. The radome on the 737 aircraft is simulated 
by the truncated fuselage as shown in Figure 5.12. 

In order to determine the appropriate polarization (horizontal or 
vertical) of the radiated field as well as the location of the antenna, 
the elevation plane patterns for a A/4 monopole, a circumferential slot 
and an axial slot mounted at Station 222 are calculated and found to be 
in very good agreement with measurements as shown in Figures 5.13(a) to 
(c), respectively. In this case, vertical polarization has better 
forward coverage than horizontal polarization because horizontal 
polarization is shorted out by the perfect conducting fuselage. 

Based on the MLS antenna coverage requirement, the radiation 
pattern of a single vertical polarization monopole is not quite 
adequate. Thus a tail mounted monopole was proposed to be mounted at 
Station 950 on the bottom of the Boeing 737 aircraft fuselage for missed 
approach purposes. This location is chosen because the slope of the 
fuselage surface at Station 950 is similar to that at Station 250 on the 
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Figure 5.12. Computer simulated model for a X/4 monopole mounted at 

Station 222 on the bottom of the fuselage of a Boeing 737 
ai rcraft. 
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(a) X/4 monopole case. 


.13. Elevation plane pattern of an antenna mounted at 

Station 222 on the bottom of a Boeing 737 aircraft. 


191 


Figure 


CALCULATED 



BOTTOM 


(SCALE: EACH OIVISION-IOOB) 


(b) Circumferential slot case. 


5.13. (Continued). 
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(c) Axial slot case. 


5.13. (Continued). 
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top of the fuselage such that the space coverage of the two antennas 
tends to have a better match. The computer model as shown in Figure 
5.14 is composed of a composite ellipsoid with four flat plates attached 
to the fuselage. Two flat plates are used to simulate the wings, and 
the other two simulate the hori zonal stabilizers. The engines are 
neglected in the present model for simplicity. The elevation plane 
radiation pattern is computed and found to be in good agreement with the 
measurements as shown in Figure 5.15. 

2. KC-135 (Boeing 707) Aircraft 

The radiation patterns for the antennas mounted on the KC-135 
aircraft are analyzed in this section. Precision pattern measurements 
(elevation and roll plane patterns) using a 1/25 scale model were taken 
at NASA (Hampton, Virginia) and will be used to verify our numerical 
solution. 

The line drawings for the aircraft are shown in Figure 5.16. The 
computer simulated models of the KC-135 for antennas mounted forward of 
the wings and over wings are shown in Figures 5.17 and 5.18, 
respectively. The elevation plane patterns for a short monopole are 
illustrated in Figure 5.19. The patterns for a circumferential KA-band 
waveguide are illustrated in Figure 5.20. Finally, the patterns for 
an axial KA-band waveguide are shown in Figures 5.21. Next, the roll 
and azimuth plane patterns for a short monopole, a circumferential 
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antenna location 


Figure 5.14. Computer simulated model for a x/4 monopole mounted at 
Station 550 on the bottom of the fuselage of a Boeing 
737 aircraft. 
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Figure 5.15. Elevation plane pattern of a X/4 monopole mounted at 
Station 950 on the bottom of the fuselage of a Boeing 
737 aircraft. 
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(a) Antenna mounted forward of wings 
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(b) Antenna mounted over wings 


Figure 5.19. Elevation plane pattern for a x/4 monopole mounted on a 
KC-135 aircraft. 
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(b) Antenna mounted over wings 

.20. Elevation plane pattern for a circumferential KA-band 
waveguide mounted on a KC-135 aircraft. 
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(a) Antenna mounted forward of wings 
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(b) Antenna mounted over wings 


>.21. Elevation plane pattern for an axial KA-band 
waveguide mounted on a KC-135 aircraft. 
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ka-band waveguide, and an axial Ka-band waveguide mounted forward and 
over the wings are shown in Figures 5.22 to 5.27. 

The computed results are found to be in very good agreement with 
the measurements in the elevation and roll planes. The comparisons of 
the azimuth plane radiation patterns are not shown here, since the 
measured data was not available. 

To further demonstrate the versatility of this solution, the 
radiation patterns for the Lindberg crossed-slot antenna [46] mounted at 
Station 470 along the top center-line of a KC-135 aircraft is 
analyzed. The Lindberg antenna (crossed slots, 90° phase difference 
between two slots) is a UHF antenna designed for use in a 
satellite-to-aircraft communications link. Using the computer model as 
shown in Figure 5.28, the radiation patterns were computed for a right 
circularly polarized Lindberg antenna. Various calculated patterns 
along with the measured results as taken form reference [46] are 
presented in Figures 5.29 to Figures 5.32 and, again, good agreement is 
obtained. The gain level in each case is adjusted to compare with 
measurements. The Eg pattern corresponds to the vertical component, 

E^ to the horizontal component and E C p pattern to the circularly 
polarized field. Note that all patterns are computed at a frequency of 
6.25 GHz and .039X x .78X slots are considered in this calculation. 
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(b) Antenna mounted over wings 


Figure 5.22. Roll plane pattern for a X/4 monopole mounted on a 
KC-135 aircraft. 
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(b) Antenna mounted over wings 


Figure 5.23. Roll plane pattern for a KA-band circumferential 
waveguide mounted on a KC-135 aircraft. 
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(b) Antenna mounted over wings 


Figure 5.24. Roll plane pattern for a KA-band axial waveguide 
mounted on a KC-135 aircraft. 
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(b) Antenna mounted over wings 


Figure 5.26. Azimuth plane pattern for a KA-band circumferential 
waveguide mounted on a KC-135 aircraft. 
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(a) Antenna mounted forward of wings 
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(b) Antenna mounted over wings 


Figure 5.27. Azimuth plane pattern for a KA-band axial waveguide 
mounted on a KC-135 aircraft. 
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(b) FRONT VIEW 



Figure 5.28. Computer simulated model for a Lindberg antenna mounted 
on a KC-135 aircraft. 
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Figure 5.29. Elevation plane pattern for a Lindberg antenna 
mounted on a KC-135 aircraft. 
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Figure 5.30. Roll plane pattern for a Lindberg antenna mounted on 
a KC-135 aircraft. 


213 


TOP 



RIGHT 

WING 


BOTTOM 

(SCRLE: ERCH DIVISION* 1 GOB) 


( c ) E cp 


Figure 5.30. (Continued). 
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Figure 5.31. Azimuth plane pattern for a Lindberg antenna mounted on a 
KC-135 aircraft. 
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(b) Eg 

Figure 5.32. Azimuthal conical pattern (9 p a 45°) for a Lindberg antenna 
mounted on a KC-135 aircraft. 
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C. PRIVATE AIRCRAFT SIMULATIONS 


This section presents a study of the effects on the radiation 
patterns associated with modeling private aircraft. 

Measured roll plane patterns were provided by NASA (Hampton, 
Virginia) [37,38,39] and are used throughout the following discussion. 

1. Cessna 402B 

Let us consider the problem of modeling the Cessna 402B illustrated 
in Figure 5.33. Our goal now is to investigate how to model the roll 
plane of this aircraft using a composite ellipsoid and flat plates to 
best match the experimental results. Experimental results have been 
obtained from NASA (Hampton, Virginia), using the 1/7 scale model at a 
range of 50 feet; however, the model is input in full scale dimensions 
at a pattern range of 350 feet. 

Past studies have indicated that matching the curvature of the 
fuselage near the antenna is more advantageous than trying to 
approximate the overall fuselage shape. With that in mind, the models 
illustrated in Figure 5.34 are considered first. Here the height and 
width of the ellipsoid match the corresponding dimensions of the actual 
fuselage top portion near the antenna, and wings are simulated as flat 
plates. Note that the prolate spheroid model of the old Aircraft Code 
[24] is also included in the figure for comparison. The resulting roll 
plane patterns are shown in Figure 5.35. The comparison shows that the 
ellipsoid model has led to a reduction in intensity of about 4 dB 
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ANTENNA LOCATION 



(a) Side view. 



(b) Front view. 


Figure 5.33. Cessna 402B. 
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(a) TOP VIEW 
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( ELLIPSOID MODEL ) ( PROLATE SPHEROID MODEL ) 

(b) FRONT VIEW 


Figure 5.34. Model for Cessna 402B with wings only. 
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throughout most of the lower half of the pattern, which is a definite 
improvement; however, the large lobes at ±105° and some ripple 
throughout the pattern are still- missing. 

A better model is obtained by consideration of the fuel tanks 
mounted on the wing tips. The actual shape of the tanks is curved, but 
of course they must be approximated by flat plates. A model with each 
fuel tank simulated as a single plate is shown in Figure 5.36, with the 
corresponding pattern shown in Figure 5.37. It is apparent that this 
is a major improvement in the lower half of the pattern, as the ±105° 
lobes are not present. A problem, however, is the excessive signal 
intensity in the 0° to ±20° range due to strong direct reflections off 
the fuel tanks. 

To alleviate this problem, let us use two plates to simulate each 
fuel tank. This will decrease the amount of direct diffraction by 
better simulating the curved surface of the tank. Such a model is 
illustrated in Figure 5.38, with the corresponding pattern shown in 
Figure 5.39. A definite improvement is shown in the figure, but the 
intensity in the ±110° to ±140° range is still about 6 or 8 dB high. 
However, this pattern is within reasonable engineering agreement with 
the experimental result. 

For greater accuracy, the effects of the engine housings can be 
considered. A model with engine housings and fuel tanks is shown in 
Figure 5.40, and the pattern resulting from this simulation model is 
shown in Figure 5.41. It is apparent from this result that the 
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Figure 5.37. Pattern for Cessna 402B model shown in Figure 5.36. 
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Figure 5.38. Model for Cessna 402B with two plate simulation of fuel 
tanks. 
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Figure 5.39. Pattern for Cessna 402B model shown in Figure 5.38. 
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Figure 5.40. Model for Cessna 402B with engines and fuel tanks. 
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Figure 5.41. Pattern for Cessna 402B model shown in Figure 5.40. 
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pattern is greatly improved in the ±120° to ±180° range. Certainly 
the results shown on this figure are acceptable for engineering 
applications. 

2. Beechcraft Baron 


Another case which presents an interesting modeling challenge is 

tho Roorhr raft Raron uji t h fho anf nrms m/Minf A/l « A * C 
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illustrated in Figure 5.42. This model is most useful for testing our 
new code in that this has wing mounted engines and propellers. 

A model is desired which will yield a computed roll plane pattern 
in close agreement with experimental results. An initial model 
including wings and engine housings is illustrated in Figure 5.43 with 
the corresponding roll plane pattern shown in Figure 5.44. One notes 
that the computed result is well matched with the experimental result, 
but the large ripple, especially in the lower half, is missing. Perhaps 
an upgrading of the model could improve our accuracy. However, the rest 
of the pattern would seem, if not perfect, at least acceptable for any 
engineering applications. 

Next, let us consider the effect of the rotating propellers in 
front of the engines. It is necessary to check the scattering due to 
the rotation of the propellers in that they are close to the antenna. 
Four different positions (i.e., 0°, 45°, 90°, 135°) of the stationary 
propellers are chosen to simulate the rotating motion of the propellers 
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Figure 5.42. Beechcraft Baron with antenna in forward location 



antenna location 


(a) TOP VIEW 



(b) FRONT VIEW 


Figure 5.43. Beechcraft Baron model with engine housings. 
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Figure 5.44. Roll conical pattern (9 p =80°) Beechcraft Baron model 
shown in Figure 5.43. 
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as shown In Figure 5.45. For simplicity, only the left propeller is 
considered with the corresponding roll plane pattern shown in Figure 
5.46. Note that width of pattern line indicates the variation of the 
radiation pattern due to the rotation of the propellers. One notes that 
the variation margin in the 200° to 300° range is about 4 or 6 dB. This 
example shows the extended capability of our code in the pattern 

^ M * k « MA J . •* i> L A i. . X J . _ .L J ^ J. 
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3. Cessna 150 

Another interesting modeling challenge is the Cessna 150 
illustrated in Figure 5.47. Perhaps the obvious way to model the Cessna 
150 would be to use the ellipsoid as a fuselage, then put in a flat 
plate above it to simulate the wing. However, the Aircraft Code has the 
requirement that the source be mounted on the ellipsoid; whereas in this 
case, the antenna is located on the wing. An acceptable model within 
the requirements of this program can be obtained as illustrated in 
Figure 5.48, where the ellipsoid is used lengthwise as the wing. Note 
that the nose and tail are modeled as flat plates. The resulting 
elevation plane pattern for this model is shown in Figure 5.49. 

Although the magnitude of ripple is not quite perfect, it is of the 
correct spacial frequency, and the general shape of the pattern is good. 

A slight disadvantage to this particular use of the ellipsoid is 
due to the sharp curvature near the antenna, in that the pattern is 
quite sensitive to small changes in source location. 
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Figure 5.45. Beechcraft Baron model with rotating properl lers on one 
side. 
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.46. Roll plane pattern (e p =80°) for Beechcraft Baron model 
shown in Figure 5.45/ 
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(b) Front view. 


5.47. Cessna 150. 
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(a) Top view. 



(b) Side view. 

Figure 5.48. Cessna 150 model. Dashed lines are not part of the 
computer simulation. 
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Figure 5.49. Cessna 150 elevation plane pattern for model shown in 
Figure 5.48. 
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D. MILITARY AIRCRAFT AND MISSILES 


This section presents the study of the radiation patterns for the 
antennas mounted on the fuselage of various military aircraft (F-16, 

F-4, A-10 and C-141) and missiles. The ability of this solution to 
analyze very complex structures is verified in terms of military 
aircraft. Various scattering mechanisms are also studied individually 
using the F-16 in order to show the significance of the variouserms. 

1. F-16 Fighter Aircraft 

The most interesting model among military aircraft studied to date 
is the F-16 fighter aircraft in that extensive experimental results [47] 
are available. The measured data was obtained by General Dynamics using 
a quarter scale model of the F-16. The line drawings for a F-16 fighter 
aircraft are shown in Figure 5.50. 

The TACAN antenna mounted on the top of the F-16 fighter fuselage 
as shown in Figure 5.50 is studied here and operated at a frequency of 
0.96 GHz. As shown in Figure 5.51, a composite ellipsoid (21.5" x 23" x 
400" x 250") is chosen to simulate the fuselage of the aircraft as 
accurately as possible, especially near the antenna location, and the 
other appendages (i.e., wings, stabilizers, tails, etc.) of the aircraft 
are simulated using 12 flat plates. Note that the radome of the F-16 
fighter is simulated as a truncated fuselage. 
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The various UTD terms are demonstrated for the roll plane radiation 
pattern in Figure 5.52 and the total field is compared with the 
experimental results in Figure 5.53. Note that the pattern of the 
prolate spheroid model [24] is also included in the figure to show the 
improvement of the new code. As can be seen from the comparison, the 
result of the ellipsoid model shows very good agreement with the maximum 
discrepancy being less than 2 dB. Note that each pattern is normalized 
to the same level such that one can compare the significance of each IJTD 
term. Since the individual UTD terms can be analyzed separately using 
this code, one can ascertain the significant features of a complex 
target. The elevation plane pattern is shown in Figures 5.54 and the 
calculated result compares quite favorably with the measurement. 

To show the complete volumetric radiation patterns, the various 
azimuthal conical patterns (0p = 10°, 20°, 30°, 40°, 45°, 50°, 55°, 

60°, 65°, 70°, 75°, 80°, 85°, 90°, 95°, 100°, 105°, 110°, 115°, 120°, 
see Figure 5.2(c)) are calculated as shown in Figures 5.55 through 5.74. 
In each case, both the principal and cross polarizations are considered. 
The calculated results compare very favorably with the measurements in 
each case. It is noted that the cockpit section simulation is not 
complete in our model; as a result, one can not expect good agreement 
between the calculated and measured results in the nose region in that 
the cockpit would be along the direct radiation path for such radiation 
directions. In addition, the ripple above the aircraft in the elevation 
pattern is most likely created by the cockpit which is not simulated in 
the analysis. 
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Figure 5.52. Various UTD terms for the roll plane pattern of a x/4 
monopole mounted on top of an F-16 fighter aircraft. 
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Figure 5.52. (Continued). 
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Figure 5.52. (Continued). 
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Figure 5.53. Total field for roll plane pattern of a X/4 monopole 
mounted on top of an F-16 fighter aircraft. 
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Figure 5.54. Elevation plane pattern of a x/4 monopole 
mounted on top of an F— 16 fighter aircraft. 
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(b) E* 

Figure 5.55. Azimuthal conical pattern (9p s 10°) of a X/4 monopole 
mounted on top of an F-16 fighter aircraft. 

(See Figure 5.2(c)). 
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Figure 5.56. Azimuthal conical pattern (Qp * 20°) of a X/4 monopole 
mounted on top of an F-16 fighter aircraft. 
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Figure 5.57. Azimuthal conical pattern (6p ■ 30°) of a X/4 monopole 
mounted on top of an F-16 fighter aircraft. 


252 


MSE 



(a) E e 


NOSE 



tnil 


(b> E* 

Figure 5.58. Azimuthal conical pattern (9 p = 40°) of a X/4 monopole 
mounted on top of an F-16 fighter aircraft. 
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Figure 5.59. Azimuthal conical pattern (9 p = 45°) of a X/4 monopole 
mounted on top of an F-16 fighter aircraft. 
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Figure 5.60. Azimuthal conical pattern (0p * 50°) of a X/4 monopole 
mounted on top of an F-16 fighter aircraft. 
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Figure 5.61. 


Azimuthal conical pattern (9p * 55°) of a A/4 monopole 
mounted on top of an F-16 fighter aircraft. 
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Figure 5.62. Azimuthal conical pattern (0 p * 60°) of a X/4 monopole 
mounted on top of an F-16 fighter aircraft. 
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Figure 5.63. Azimuthal conical pattern ( 9p = 65°) of a X/4 monopole 
mounted on top of an F-16 fighter aircraft. 
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Figure 5.64. Azimuthal conical pattern (9p = 70°) of a X/4 monopole 
mounted on top of an F-16 fighter aircraft. 
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Figure 5.65. Azimuthal conical pattern (9 p * 75°) of a X/4 monopole 
mounted on top of an F-16 fighter aircraft. 
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.66. Azimuthal conical pattern (9p = 80°) of a A/4 monopole 
mounted on top of an F-16 fighter aircraft. 
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Figure 5.67. Azimuthal conical pattern (Qp » 85°) of a X/4 monopole 
mounted on top of an F-16 fighter aircraft. 
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Figure 5.68. Azimuthal conical pattern (0p » 90°) of a X/4 monopole 
mounted on top of an F-16 fighter aircraft. 
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Figure 5.69. Azimuthal conical pattern (0 p * 95°) of a X/4 monopole 
mounted on top of an F-16 fighter aircraft. 
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Figure 5.70. Azimuthal conical pattern (9n a 100°) of a X/4 monopole 
mounted on top of an F-16 fighter aircraft. 
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Figure 5.71. Azimuthal conical pattern (9 p * 105°) of a X/4 monopole 
mounted on top of an F-16 fighter aircraft. 
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Figure 5.72. Azimuthal conical pattern (9p * 110°) of a X/4 monopole 
mounted on top of an F-16 fighter aircraft. 
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Figure 5.73. Azimuthal conical pattern (0p a 115°) of a X/4 monopole 
mounted on top of an F-16 fighter aircraft. 


268 


ROSE 



TRIL 


(b) E + 


Figure 5.74. Azimuthal conical pattern (9p * 120°) of a X/4 monopole 
mounted on top of an F-16 fighter aircraft. 
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2. F-4 Fighter Aircraft 


This example is illustrative in that the antenna is mounted on the 
belly of an F-4 aircraft which is loaded with armament. The F-4 fighter 
aircraft has a complex airframe with many scattering objects (i.e., fuel 
tanks, jet intakes, pylons, etc.) attached to the basic frame. A X/4 
monopole is used to simulate a UHF blade which is mounted on the bottom 
fuselage of an F-4 aircraft and operated at a frequency of 375 MHz. The 
measured data was obtained at the RADC Newport antenna range. 

The line drawings of the F-4 fighter aircraft are shown in Figure 

5.75. The computer model of the F-4 aircraft is illustrated in Figure 

5.76. A composite ellipsoid (5" x 20" x 300" x 200") is chosen to 
simulate the bottom surface of the aircraft fuselage; whereas, the other 
appendages are simulated using 11 flat plates. It is noted that the 
actual shape of the fuel tank is curved and each fuel tank is simulated 
as a single plate in our model for simplicity. 

The azimuthal conical pattern (9p * 105°) is compared with a 
measurement result in Figure 5.77. Although some discrepancy does 
exist, the general shape of the pattern is in good agreement. 


270 




Figure 5.76. Computer simulated model of an F-4 Phantom fighter 
aircraft. 


Ill 


< 


CALCULATED 

MEASURED 


NOSE 



TAIL 

(SCALE: EACH DIVISION. 10DB) 


Figure 5.77. Azimuthal conical pattern ( 0 p s 105°) of a X/4 monopole 
mounted on the belly of an F-4 fighter aircraft. 
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3. A-10 Aircraft 


In this case, a four monopole array is mounted on the belly of an 
A-10 aircraft. The line drawings are shown in Figure 5.78. The 
experimental data was also obtained at the RADC Newport site. 

The computer simulation model is illustrated in Figure 5.79. In 
that the four monoples are spaced a half wavelength apart, the mutual 
coupling between the antennas is significant and cannot be neglected in 
the pattern calculations. This can be solved using a thin wire moment 
method code [49,50] to treat the four closely spaced loaded-dipoles (50 
ohms at center of each dipole) in order to obtain the excitation of each 
monopole even though only one element was excited. 

The azimuthal conical pattern (0p = 105°) is compared with 
measurement in Figure 5.80. The result reveals good engineering 
agreement. 

4. C-141 Aircraft 

Let us consider a monopole mounted on the top of a C-141 aircraft. 
The line drawings for a C-141 aircraft are shown in Figure 5.81. The 
experimental work [40] was performed at General Dynamics (San Diego, 
California) using a 1/10 scale model of a C-141 aircraft. 

Figure 5.82 illustrates the computer model which is used to 
simulate a C-141 aircraft. A composite ellipsoid (7. 37' x 8.37' x 90' x 
46.05') is chosen to simulate the fuselage surface of the actual 
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(a) FRONT VIEW 



(b) BOTTOM V I EW 


Figure 5.79. Computer simulated model of an A-10 aircraft. 
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Figure 5.80. Azimuthal conical pattern (8 p * 105°) of four monopoles 
mounted on the belly of an A-10 aircraft. 
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(c) TOP VIEW 


Figure 5.81. C-141 aircraft. 
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(a) SIDE VIEW 



(b) FRONT VIEW 



(c) TOP VIEW 


Figure 5.82. Computer simulated model of a C-141 aircraft. 
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aircraft as accurately as possible, and finite flat plates are used to 
simulate the wings, vertical stabilizer and T-tail as shown in Figure 
5.82. 

The radiation patterns for a X/4 monopole mounted on the fuselage 
of a C-141 aircraft are calculated using the model just described. In 
order to verify the validity of the analytic solution, various azimuthal 
conical patterns (10°, 20°, 30°, 40°, 50°, 60°, 70°, 80°, 90°, 100°, see 
Figure 5.2(c)) and the elevation plane pattern are computed and compared 
with measured results as shown in Figures 5.83 and 5.84, respectively. 
The calculated results compare very favorably with the scale model 
measurements. 

5. Missile 

The missile antenna problem considered next is interesting in that 
the antenna is mounted among four large ram jets which significantly 
impact on the antenna location problem. 

The front view of the missile and the computer simulation model are 
shown in Figure 5.85. Note that the axial slot antenna is mounted 
between the ram jets. The fuselage is simulated as a composite 
ellipsoid and only two side walls of the intakes are simulated with two 
flat plates. It is noted that a fuselage blockage plate is added in the 
computer simulation model and is shown as a dashed line in Figure 
5.85(b). This feature is added in the code to prevent the ray 
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Figure 5.83. Azimuthal conical patterns of a X/4 monopole mounted on 
a C-141 aircraft. 
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Figure 5.83. (Continued). 
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Figure 5.83. (Continued). 
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Figure 5.83. (Continued). 
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Figure 5.84. Elevation plane pattern of a X/4 monopole mounted on 
top of a C— 141 aircraft. 
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(a) MISSILE MODEL 



(b) COMPUTER SIMULATED MODEL 


Figure 5.85. Missile model for an axial slot mounted between two 
ram jets. 
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contribution passing through the fuselage itself in the field 
calculation. The roll and elevation plane patterns of this model are 
compared with measured results in Figures 5.86 and 5.87, respectively. 
The calculated results are in good agreement with the measurement. 

In order to improve the pattern performance, the axial slot antenna 
was mounted on the ram jet instead of the fuselage as shown in Figure 
5.88. In this case, the air duct is simulated as a composite ellipsoid 
and the other structures (including the fuselage) are simulated by 
multiple flat plates. The roll and elevation plane patterns are 
compared with measurement data in Figures 5.89 and 5.90, respectively. 
The results reveal very good agreement between the theoretical 
prediction and actual measurements. 

E. SPACE SHUTTLE 

In the design of antennas for spacecraft, re-entry vehicles, and 
high performance aircraft, the major area of difficulty is achieving 
near omni-directional pattern coverage with flush mounted elements. The 
problem is compounded by the unavailability of several prime locations 
on the frame, such as the top of the mid fuselage, all leading-edge 
surfaces, the nose, and the vertical stabilizer. These difficulties 
lead to the Orbiter antennas being specially developed and qualified in 
order to meet the performance requirements with a high degree of 
reliability. In addition, the radiation coverage requirements of some 
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Figure 5.86. Roll plane pattern for an axial slot mounted between two 
ram jets. 
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Figure 5.87. Elevation plane pattern for an axial slot mounted between 
two ram jets. 
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(b) COMPUTER SIMULATED MODEL 


Figure 5.88. Missile model for an axial slot mounted on a ram jet. 
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.89. Roll plane pattern for an axial slot mounted on a 
ram jet. 
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Figure 5.90. Elevation plane pattern for an axial slot mounted on a 
ram jet. 
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of the Orbiter antennas are complicated by the multiple requirements for 
operation during launch, from earth orbit to both ground station and 
relay satellites, and with the landing area during atmospheric flight. 

The major factor influencing the decision as to what type of 
antenna design to utilize for each of the Shuttle systems is the 
determination of the mounting location on the Shuttle vehicle. The 
choice of a mounting location is influenced by such factors as available 
mounting structure, TPS (thermal protection system) covering material, 
required angular coverage, and the effects of nearby blocking structure 
during various phases of the Shuttle mission timeline. The task of 
determining antenna locations is normally accomplished utilizing a 
full-scale mockup of the subject vehicle on an antenna range. However, 
due to the large size of the Shuttle, such an approach is not feasible. 
Therefore, the need for an efficient analytic solution is quite 
apparent. 

The performance of the S-Band Quad Antennas for PM communications 
is investigated here using this analytic solution. The basic antenna 
configuration [48] consists of a set of four circularly polarized 
antennas flush-mounted with 90° spacing around the roll axis of the 
Orbiter forward fuselage as shown in Figure 5.91. 

The computer simulation model for the top mounted antenna is 
illustrated in Figure 5.92. A composite ellipsoid (75" x 104" x 1500" x 
130") is chosen to simulate the top portion of the fuselage surface, 
and the other structures such as side walls, nose, wings and vertical 
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Figure 5.91. S-band quad antenna locations on a Space Shuttle. 


295 




(a) 3-DIMENSIONAL VIEW 




Figure 5.92. Computer simulated model for a crossed-slot antenna 
mounted on top of a Space Shuttle Orbiter. 
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stabilizer are simulated by 23 flat plates as shown in Figure 5.92. The 
three principal plane radiation patterns for the top mounted 
crossed-slot antenna are computed as shown in Figures 5.93 through 5.95. 

The blocking effects of the payload bay doors during various phases 
of the Shuttle mission timeline are analyzed by calculating the roll 
conical (e p =45°) patterns with the doors open and closed. For the top 
mounted antenna, the heat radiators which actually cause the blocking 
and reflection of the incident ray are simulated as three flat plates on 
each side as shown in Figure 5.96. The roll conical (6p=45°) patterns 
for the shuttle with payload bay doors opened are compared with those 
for the closed payload doors in Figures 5.97 and 5.98. Comparing these 
patterns, one can easily see the blocking effect of the heat 
radiator in the 100° to 160° range and the large disturbance (around 20 
dB) in the pattern in 30° to 60° range which is the main beam region of 
the top mounted antenna. This is due to strong direct reflections off 
the heat radiator. 

The computer simulation model for the bottom mounted antennas is 
illustrated in Figure 5.99. Note that the simulation models are turned 
upside down for bottom mounted antennas. For this situation, the 
outer covers of the payload bay doors cause the blockage and reflection, 
and are simulated as two flat plates on each side as shown in Figure 
5.100. The roll conical patterns (0p=45°) for the shuttle with 
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Figure 5.93. Roll plane radiation patterns for a crossed-slot antenna 
mounted on top of a Space Shuttle Orbiter. 
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.94. Elevation plane patterns for a crossed-slot antenna 
mounted on top of a Space Shuttle Orbiter. 
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Figure 5.95. Azimuth plane patterns for a crossed-slot antenna 
mounted on top of a Space Shuttle Orbiter. 
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(c) TOP VIEW 


Figure 5.96. Computer simulated model for a crossed-slot antenna 
mounted on top of a Space Shuttle Orbiter when the 
payload bay doors are open. 
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Figure 5.97. Roll conical patterns (9g s 45°) for a crossed-slot 
antenna mounted on top of a Space Shuttle Orbiter 
when the payload bay doors are closed. 
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Figure 5.98. Roll conical patterns (00*45°) for a crossed-slot 
antenna mounted on top or a Space Shuttle Orbiter 
when the payload bay doors are open . 
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Figure 5.99. 


Computer simulated model for a crossed-slot antenna 
mounted on the bottom of a Space Shuttle Orbiter. 




la] SIDE VIEW 


Figure 



(b) FRONT VIEW 



.100. Computer simulated model for a crossed-slot antenna 
mounted on the bottom of a Space Shuttle Orbiter when 
the payload bay doors are open. 


payload bay doors opened are compared with those closed in Figures 5.101 
and 5.102, respectively. The results also indicate severe disturbance 
of the pattern due to the strong, reflections off the payload bay doors. 

The above results indicate why one can lose use of his 
communication channel during those times when the shuttle payload doors 
are open. This illustrates how one can quickly analyze the performance 
of antennas in a complex changing environment. 
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Figure 5.101. Roll conical patterns (9 p =45°) for a crossed-slot antenna 
mounted on the bottom of a Space Shuttle Orbiter when 
the payload bay doors are closed. 
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Figure 5.102. Roll conical platterns (0n=45°) for a crossed-slot 
antenna mounted on the bottom of a Space Shuttle 
Orbiter when the payload bay doors are open . 
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CHAPTER VI 


SUMMARY AND CONCLUSIONS 

The object of this research has been to develop an efficient and 
accurate analytic solution for the high frequency radiation patterns of 
on-aircraft antennas mounted on the fuselage. 

In most cases, scale model measurements have been used to design, 
and locate airborne antennas. However, this approach is both expensive 
and time consuming. Thus, the need for an efficient and reliable 
analytic solution of on-aircraft antenna patterns is quite apparent. 

The approach applied for this solution is the Uniform Geometrical Theory 
of Diffraciton (UTD). It is a high frequency technique with the only 
limitation being that the source and various scattering centers be 
separated on the order of a wavelength or more. In some cases, this 
requirement can be relaxed to approximately a quarter-wavelength. 

The basic approach applied here is to break up the aircraft into 
its simplest structural form so that the study is applicable to 
general -type aircraft. The fuselage is modeled as a perfectly 
conducting composite ellipsoid; whereas, the other appendages (such as 
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wings, horizontal and vertical stabilizers, nose, fuel tanks and 
engines, etc.) are modelled as perfectly conducting flat plates that can 
be attached to the fuselage or to each other. The composite ellipsoid 
fuselage model is necessary to successfully simulate the wide variety of 
real world fuselage shapes. In fact, the fuselage has a dominant effect 
on the resulting radiation pattern in that the significant energy region 
is confined to the vicinity of the antenna. 

An efficient numerical technique, which uses elliptic cylinder 
and elliptic cone perturbation methods, for the computation of the 
geodesic paths on an ellipsoid surface was introduced and verified by 
the comparisons with exact solutions in Chapter III. The results 
confirm that for a given radiation direction in the shadow region, the 
geodesic path (i.e., the actual ray path) and the final diffraction 
point on the ellipsoid surface can be efficiently solved using the 
perturbation methods. Based on the fuselage simulation techniques 
discussed in Chapter III, the algorithms for the actual computation of 
each individual UTD term such as the source field, reflected field, 
diffracted field, corner diffracted field, and higher order terms have 
been discussed in Chapter IV. 

Using the fuselage simulation technique and the field calculation 
algorithms, this new Aircraft Code has been developed to compute and 
superimpose the various UTD terms for near or far field calculations. 
This newly developed analytic solution provides an accurate, efficient 
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and useful means for determination of more optimum antenna designs and 
locations based on their pattern performance. The capability of the new 
solution was demonstated for the- wide variety of examples in Chapter V. 

The optimum antenna locations and proper polarization have been 
found for the microwave landing system on a Boeing 737 aircraft. The 
scattering effects of the wings, tail, engine housing, fuel tanks, and 
rotating propellers on the radiation patterns were analyzed for the 
various antennas mounted on a KC-135, a Cessna 402B and a Beechcraft 
Baron. Note that the rotating motion of the propellers was simulated by 
considering 4 different stationary positions of the propellers. The 
program runs a pattern in approximatley 4 minutes for a commercial 
aircraft model on a VAX 11/780 computer. 

To show the applicability of this solution for an antenna mounted 
on the wing of a Cessna 150 aircraft, the ellipsoid was used lengthwise 
as the wing, and the fuselage was modelled by flat plates. 

To further demonstrate the versatility of this solution, various 
azimuthal radiation patterns for a monopole mounted on an F-16 fighter 
and C-141 aircraft were obtained. The scattering effects of complex 
stores mounted on military aircraft were also investigated in terms of 
an F-4 fighter and A-10 aircraft. Missile antennas mounted among four 
large ram jets were also simulated for both an antenna mounted between 
as well as on the ram jets. 
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The validity of this solution was also verified by numerous 
comparisons between the calculated and measured data. In most cases, 
the results of this solution revealed an improvement over the previous 
solutions [24,40] and show very good agreement between the theoretical 
predictions and various scale model measurements. 

In a more recent example, the performance of the S-Band Quad 
antennas mounted on the Space Shuttle Orbiter was anal y zsd . The 
scattering and blocking effects of the open payload bay doors was 
analyzed by computing the roll conical patterns (e p =45°). The results 
revealed large disturbances (more than 20 dB) in the pattern in the main 
beam regions for both the top and bottom mounted antennas. This is due 
to the strong reflections off the payload bay doors and heat radiators. 
Using this approach, one can efficiently analyze the pattern performance 
of antennas radiating in a complex environment such as an aircraft. 

With such a tool, an antenna designer can consider numerous alternatives 
in order to optimize his design. 
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APPENDIX A 


EXACT GEODESIC PATHS FOR AN ELLIPSOID 

A numerical method employing calculus of variations to calculate 
the geodesic paths on an ellipsoid C36] is presented in this appendix. 
This solution is used to check the validity of the perturbation solution 
which was introduced to efficiently calculate the geodesic paths on an 
ellipsoid in Section III-C. 

Using rectangular coordinates, an ellipsoid can be described as 
x 2 y 2 z 2 

72 + b 2 + c 2 = 1 (A,1) 

where, without loss of generality, it is assumed that c>b>a>0. The key 
to the derivation of the geodesic path solution of an ellipsoid is to 

find a coordinate system which is orthogonal on the ellipsoidal surface. 


313 



Consider the following three equations 


x 2 y 2 z 2 

m + p^ + ^ sl * a2 ' > 5 


(A.2) 



b 2 > n > a 2 


(A.3) 



y 2 z 2 
b 2 -r + c 2 -t 


= 1 


c 2 > n > b 2 


(A. 4) 


which are of an ellipsoid, a hyperboloid of one sheet, and a hyperbolid 
of two sheets [35], respectively, all confocal with the ellipsoid of 
Equation (A.2). The variables u 1 = c, u 2 = n, u 3 = t are called 
ellipsoidal coordinates. The transformation to the rectangular 
coordinates is obtained by solving Equations (A.2), (A.3), and (A. 4) 
simultaneously for x, y and z, such that 


— (a 2 -0(a 2 -n)(a 2 -T)~ 
(c 2 -a 2 )(b 2 -a 2 ) 


1/2 


(A. 5) 


= + 


-(b 2 -0(b 2 -n)(b 2 - T )- 

(c 2 -b 2 )(a 2 -b 2 ) 


1/2 


(A. 6) 


z = ± 


“(c 2 -C)(c 2 -n)(c 2 -x)— 

(a 2 -c 2 )(b 2 -c 2 j 


1/2 


(A. 7) 
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In terms of the ellipsoidal coordinates, the displacement vector 


dr can be written as 


3r 3? 3r 

dr » du 1 + -^2 du2 + du3 s V ul + h duZ + ®3 dlj3 


Then the length of a line element, denoted by ds, is 


3 3 

(ds)^ = dr • dr * £ I dlj1 du ^ 

i=l j=l 1 J 


* l I 9jj du 1 du j 
i=l j=l 


(A. 


(A. 


where 


3x 3 x 3y 3y 3z 3z 

9 ij = a i a j 3U 1 3uJ + 3U 1 3u3 + 3U 1 3uJ 

It can be shown, using Equations (A. 5), (A. 6) and (A. 7), that 
g i j ■ o, if i * j 

1 (T-O(n-e) 

g U = A (a 2. c)(b 2. 5 )( c 2_ c) 

1 U-n)(T-n) 

g 22 = 4 (a 2 -n)(b 2 -n)(c 2 -n) 

and 

1 (n-x)(g-t) 

g 33 a 4 (a 2 -T)(b 2 -T)(c 2 -t) 


(A. 

(A 

(A 

(A 

(A 



By substituting Equations (A. 11) through (A. 14) into Equation (A. 9), one 
obtains 


(ds) 2 = 


l 

4 


(t-0(n-0 ' 
(a 2 -S)(b 2 -S)(c 2 -0 


2 (5-n)(x-n) 

(d5) + ( a 2_ n ) (b 2 -n) (c 2 -n) 


(dn) 2 


(n-x)(5-x) ? 

+ (a 2 -x)(b 2 -x)(c 2 -x) (dT) 


(A. 15) 


For the ellipsoid of Equation (A.l), S = d? * 0. Then Equation (A. 15) 
becomes 


(ds) 2 


4 (x-n) 


n(dn) 2 

(n-a 2 )(Ti-b 2 )(n-c 2 ) 


T(d t) 2 

( t-a 2 ) ( t-b 2 ) ( t-c 2 ) 


(A. 16) 


or 


1/2 


1 n s 1/2 

s - 2 /(x- n) 


nn 


•2 


(n-a 2 )(n-b 2 )(ri-c 2 ) ‘ (x-a 2 )(x-b 2 )(x-c 2 ) 


dt 


(A. 17) 


where n 1 s dn/dx and dx is assumed to be positive. Let h(n,n',x) denote 
the integrand of Equation (A. 17), i.e.. 
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and from Equation (A. 19), one obtains 

d 3h 2 dh 3h 3h 2 
dx 3n' ^ dx 3n' + 3n 


By substituting Equations (A, 19) and (A. 20) into Equation (A. 21) and 
rearranging the various terms, one obtains 



d[ T|* ( T-n)] 1 3h 


Id 3h 

"dr h Tn 7 + n ‘ ^ T_n ^ h dr TV 


1 dh 3h 

n'( T -n) hi dTT7 + n' - 0 


which can be simplified to 


A 

St 

dx 


1 3h 

n ' (T - , ' ) h 77 * 


(n-b z ) 


= 0 


(A. 22) 


(A. 23) 


Equation (A. 19) becomes a complete differential. The derivation of the 
geodesic path solution is straight forward from Equation (A. 23). It is 
obvious that 

1 3b 

VC T-n) -jj- *^T + n - b = -B (A. 24) 


where e is a constant. 

Combining Equations (A. 18) and (A. 24), one obtains 

- o(dn) 2 - T(dx ) 2 

(n-a 2 )(n-b 2 )(n-c 2 )(o-b 2 +8) ” (x-a 2 )( t-b 2 ) ( t-c 2 ) ( t-b 2 +B) * (A. 25) 

Since c 2 >t>b 2 >n>a 2 , it is obvious that T-b 2 +e>0>n-b 2 +8 or b 2 -n>S>b 2 -x. 
Thus, Equation (A. 25) is the geodesic path solution. However, it is 
more convenient to make the following changes of variables. 
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n 


a 2 si n 2 4> + b 2 cos 2 4> , 


(A. 26) 


and 

t * b 2 cos 2 4> + c 2 sin 2 *p . (A. 27) 

In terms of $ and f, the geodesic path solution in Equation (A. 25) 
can be rewritten as 

[b 2 - (b 2 -a 2 ) sin 2 »] (d») 2 

[(c 2 -b 2 ) + (b 2 -a 2 ) sin 2 <J>] [-6 + (b 2 -a 2 ) sin 2 <|>] 

[b 2 + (c 2 -b 2 ) sin 2 tl (d4») 2 

3 [(b 2 -a 2 ) + (c 2 -b 2 ) sin 2 i|»] [6 + (c 2 -b 2 ) sin 2 *] * (A. 28) 


and in integral form as 


1/2 


♦ 

/ 

*1 


[b 2 - (b 2 -a 2 ) sin 2 <f>] 


[(c 2 -b 2 ) + (b 2 -a 2 ) sin 2 <|>] [-6 + (b 2 -a 2 ) sin 2 <fr] 


|d<fr| 


1/2 


/ 

♦l 


b 2 + (c 2 -b 2 ) sin 2 t 


[(b 2 -a 2 ) + (c 2 -b 2 ) sin 2 *] [6 + (c 2 -b 2 ) sin 2 *] 


|d*| 


(A. 29) 


where b 2 -a 2 >S>b 2 -c 2 . Note that the absolute values of di> and d\p are 
used in Equation (A. 29). 
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To define the ranges of $ and *, it is necessary to go back to 
Equations (A. 5), (A. 6) and (A. 7). With £=0 and in terms of 4> and 
they can be rewritten as 

d 1/2 

X - ( c 2. a 2)i/2 cos * {b 2 zosH + c 2 sin 2 * - a 2 ) (a. 30) 

y = bsi n <t> sin* (A. 31) 


and 

2 ■ ( c 2 _ a 2 )i /2 cos* (c 2 -a 2 s1n* + - b 2 cos 2 *) 1/2 . (a. 32) 

If the goedesic path. Equation (A. 29), crosses the curve **0 or ir, then 
3<0 and if it crosses the curve *=0 or ir, then 8>0. Thus to ensure the 
continuity of * and * along the geodesic path, the ranges of * and * are 
defined as follows: 

0 < * < it, -ir < * < 2ir for 8 < 0 


and 

-ir < * < 2ir, 0 < * < it for 8 > 0 . 

Figure A.l illustrates the * and * curves as projected onto the x-z 
plane. Equation (A, 29) is the geodesic path solution of an ellipsoid 
which is employed for obtaining more exact geodesic paths. When the 


320 



Figure A 


t 

^*0 



t p i (*o»yo.* 0 ) 


4 > * o 


->> x 


<£= CONSTANT 

^ = CONSTANT 

x 


m(£4)* 
• □ 


y a *o 


-■«(SeS)* 


.1. Projection of 4>- and ip- curves onto the xz plane of 
an ellipsoid. 
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geodesic path includes one of those four points, P^, Pg, ^3 anc * P4 
the figure, then 8 must equal zero and the integrals in Equation (A. 29) 
diverge at those points. Thus if (<t>g,'l'g) = ^i * ^ = then it 

necessary to replace (♦g*^ ^ (<l>g + A<k, w ^ ere Ai}> » 0 and Aip - 0. 

The geodesic path between the two points and ( <i>2 » ^2 ^ can 

determined from Equation (A. 29) by first determining the value of 8. 
Since there are absolute value-signs attached to the differentia! d<f> and 
dip in the geodesic path solution, it is important to know how the 
variables <j> and i|> change from to $g an< ^ ’f'j to 'kg* respectively. 


322 



REFERENCES 


[1] J.B. Keller, "Geometrical Theory of Diffraction", J. Opt. Soc. 

Am., Vol . 52, No. 2, pp. 116-130, 1962. 

[2] P.H. Pathak, N. Wang, W.D. Burnside and R.G. Kouyoumjian, "A 
Uniform GTD Solution for the Radiation from Sources on a Convex 
Surface", IEEE Trans, on Antennas and Prop., Vol, AP-29, No. 4, 
pp. 609-622, July 1981. 

[3] R.G. Kouyoumjian and P.H. Pathak, "A Uniform Geometrical Theory of 
Diffraction for an Edge in a Perfectly-Conducting Surface", Proc. 
IEEE, Vol. 62, November 1974, pp. 1448-1461. 

[4] R.G. Kouyoumjian, P.H. Pathak and W.D. Burnside, "A Uniform GTD 
for the Diffraction by Edges, Vertices and Convex Surfaces", book 
chapter in Theoretical Methods for Determining the Interaction of 
Electromagnetic Waves with Structures , edited by J.K. Skwirynski, 
Sijthoff and Noordhoff, Netherlands, 1978. 

[5] W.D. Burnside, N. Wang and E.L. Pelton, "Near-Field Pattern 
Analysis of Airborne Antennas", IEEE Trans, on Antennas and Prop., 
Vol. AP-28, pp. 318-327, May 1980. 

[6] F.A. Sikta, W.D. Burnside, T.T. Chu and L. Peters, Jr., 
"First-Order Equivalent Current and Corner Diffraction Scattering 
from Flat Plate Structures", IEEE Trans, on Antennas and Prop., 
Vol. AP-31, No. 4, pp. 584-589, July 1983. 

[7] W.D. Burnside, "Analysis of On-Aircraft Antenna Patterns", 
Technical Report 3390-1, August 1972, The Ohio State University 
ElectroScience Laboratory, Department of Electrical Engineering; 
prepared under Contract 62269-72-C-0354 for the Department of 
the Navy, Naval Air Development Center. Also a Ph.D. 

Dissertation, The Ohio State University, 1972. 

[8] W.D. Burnside, M.C. Gilreath, R.J. Marhefka and C.L. Yu, "A Study 
of KC-135 Aircraft Antenna Patterns", IEEE Trans, on Antennas and 
Prop., Vol. AP-23, pp. 3090-316, May 1975. 

[9] C.L. Yu, W.D. Burnside, and M.C. Gilreath, "Volumetric Pattern 
Analysis of Airborne Antennas", IEEE Trans, on Antennas and Prop., 
Vol. AP-26, pp. 636-641, September 1978. 


323 



[10] C.C. Huang, N. Wang, and W.D. Burnside, "The High-Frequency 
Radiation Patterns of a Spheroid-Mounted Antenna", Report 
712527-1, March 1980, The Ohio State University ElectroScience 
Laboratory, Department of Electrical Engineering; prepared under 
Contract No. N00019-80-C-00'50 for the Department of the Navy, 

Naval Air Systems Command. 

[11] W.D. Burnside and H.H. Chung, "Analysis of Airborne Antenna 
Pattern and Mutual Coupling and their Effects on Adaptive Array 
Performance", Technical Report 714258-5, June 1983, The Ohio State 
University ElectroScience Laboratory, Department of Electrical 
Engineering; prepared under Contract No. F30602-182-C-0009 for the 
Department of the Air Force, Rome Air Development Center. Also a 
Ph.D. dissertation. The Ohio State University, 1983. 

[12] W.D. Burnside and N. Wang, "Near Field Analysis of Airborne 
Antennas", Final Report 713321-4, December 1981, The Ohio State 
University ElectroScience Laboratory, Department of Electrical 
Engineering; prepared under Contract No. N00019-80-C-0593 for the 
Department of the Navy, Naval Air Systems Command. 

[13] J.J. Kim, N. Wang, and C.D. Chuang, "Geodesic Paths of an 
Ellipsoid-Mounted Antenna", Technical Report 713321-3, March 1982, 
The Ohio State University ElectroScience Laboratory, Department 

of Electrical Engineering; prepared under Contract No. 
N00019-80-PR-RJ015 for the Department of the Navy, Naval Air 
Systems Command. 

[14] J.J. Kim and W.D. Burnside, "Geodesic Paths for Side-Mounted 
Antenna on an Ellipsoid Model", Technical Report 714215-1, 

October 1982, The Ohio State University ElectroScience Laboratory, 
Department of Electrical Engineering; prepared under Contract No. 
N00019-81-C-0424 for the Department of the Navy, Naval Air Systems 
Command . 

[15] W.D. Burnside, J.J. Kim and B. Grandchamp, "Airborne Antenna 
Radiation Pattern Code User’s Manual", Technical Report 716199-2, 
December 1984, The Ohio State University ElectroScience Laboratory, 
Department of Electrical Engineering; prepared under Grant NSG 1498 
for NASA, Langley Research Center, Hampton, Virginia. 

[16] Short Course Notes on "The Modern Geometrical Theory of 
Diffraction", Vols. 1, 2 and 3, The Ohio State University 
ElectroScience Laboratory, Department of Electrical Engineering, 
1980. 

[17] M.M. Lipschutz, Differential Geometry, Schaum’s Outline Series. 
McGraw-Hill Book Co., New York, New York, 1969. 


324 



[18] W.D. Burnside, N. Wang, and E.L. Pelton, "Near Field Pattern 
Computations for Airborne Antenna", Technical Report 784685-4, 

June 1978, The Ohio State University ElectroScience Laboratory, 
Department of Electrical Engineering; prepared under Contract No. 
N00019-77-C-0299 for the Department of the Navy, Naval Air Systems 
Command. 

[19] P.H. Pathak, and R.G. Kouyoumjian, "The Dyadic Diffraction 
Coefficient for a Perfectly Conducting Wedge", Technical Report 
2183-4, June 1970, The Ohio State University ElectroScience 
Laboratory, Department of Electrical Engineering; prepared under 
Contract No. AF19 ( 628)-5929 for the Department of the Air Force, 
Air Force Cambridge Research Laboratories. 

[20] R.G. Kouyoumjian and P.H. Pathak, "The Dyadic Diffraction 
Coefficient for a Curved Edge", Technical Report 3001-3, August 
1973, The Ohio State University ElectroScience Laboratory, 
Department of Electrical Engineering; prepared under Grant No. 

NGR 36-008-144 for the National Aeronautics and Space 
Administration. 

[21] W.D. Burnside and P.H. Pathak, "A Corner Diffraction Coefficient", 
to appear. 

[22] J.J. Kim and W.D. Burnside, "Radiation Patterns of an Antenna 
Mounted on the Off-Mid Section of an Ellipsoid", Technical Report 
714215-3, June 1983, The Ohio State University ElectroScience 
Laboratory, Department of Electrical Engineering; prepared under 
Contract No. N00019-81-C-0424 for the Department of the Navy, 

Naval Air Systems Command. 

[23] J.J. Kim and W.D. Burnside, "Radiation Patterns of an Antenna 
Mounted on the Mid Section of an Ellipsoid", Technical Report 
714215-2, July 1983, The Ohio State University ElectroScience 
Laboratory, Department of Electrical Engineering; prepared under 
Contract No. N00019-81-C-0424 for the Department of the Navy, 

Naval Air Systems Command. 

[24] H.H. Chung and W.D. Burnside, "General 3D Airborne Radiation 
Pattern Code User's Manual", Technical Report 711679-10, July 
1982, The Ohio State University ElectroScience Laboratory, 
Department of Electrical Engineering; prepared under Contract No. 
F30602-79-C-0068 for the Department of the Air Force, Rome Air 
Development Center, 


325 



[25] E.L. Pelton and W.D. Burnside, "Edge Diffraction Point Analysis 
Used in Near-Field On-Aircraft Antenna Studies", Report 784685-2, 
October 1977, The Ohio State University ElectroScience Laboratory, 
Department of Electrical Engineering; prepared under Contract No. 
N00019-77-C-0299 for the Department of the Navy, Naval Air Systems 
Command. 

[26] S.W. Lee, "Differential Geometry for GTD Application", Report No. 
77-21, October 1977, University of Illinois at Urbana-Champaign, 
Department of Electrical Engineering, Electromagnetics Laboratory; 
supported by the National Science Foundation. 

[27] R.G. Kouyoumjian, Numerical and Asymptotic Techniques in 
Electromagnetics, Topics in Applied Physics, Volume 3, edited by 
R. Mittra, Springer-Verlag, New York, pp. 164-215, 1975. 

[28] R.G. Kouyoumjian, "Advanced Electromagnetics", The Ohio State 
University EE 819 Class Notes, 1981. 

[29] T.T. Chu, "First Order Uniform Geometrical Theory of Diffraction 
Analysis of the Scattering by Smooth Structures", Ph.D. 
Dissertation, The Ohio State University, Department of Electrical 
Engineering, 1983. 

[30] W.D. Burnside, "Microwave Optics (Geometrical Theory of 
Diffraction)", The Ohio State University EE 712 Class Notes, 1981. 

[31] C.E. Ryan and L. Peters, Jr., "Evaluation of Edge-Diffracted 
Fields Including Equivalent Currents for the Caustic Regions", 

IEEE Trans, on Antennas and Prop., Vol . AP-17, pp. 292-299, May 
1969; see also correction to this paper in Vol. AP-18, p. 275, 
March 1970. 

[32] R.F. Harrington, Time Harmonic Electromagnetic Fields , McGraw-Hill 

Book Co., New York, 1961. ' ~ 

[33] C.A. Balanis, Antenna Theory Analysis and Design , Harper & Row 
Publishers, New York, 1982. 

[34] R.C. Hansen, ed., Geometric Theory of Diffraction , IEEE Press, 
1981. 

[35] J.A. Stratton, Electromagnetic Theory , McGraw-Hill Book Co., 

New York, 1941. 

[36] C.D. Chuang, "The Geodesic Path Solution of an Ellipsoid", 
unpublished paper. 


326 



[37] W.D. Burnside and K.W. Burgener, "Analysis of Private Aircraft 
Antenna Patterns', Report 710964-1, September 1978, The Ohio State 
University ElectroScience Laboratory, Department of Electrical 
Engineering; prepared under Grant NSG 1498 for National 
Aeronautics and Space Administration, Langley Research Center, 
Hampton, Virginia. 

[38] K.W. Burgener and W.D. Burnside, "Analysis of Private Aircraft 
Antenna Patterns", Report 710964-2, January 1979, The OhioState 
University ElectroScience Laboratory, Department of Electrical 
Engineering; prepared under Grant NSG 1498 for National 
Aeronautics and Space Administration, Langley Research Center, 
Hampton, Virginia. 

[39] W.D. Burnside, "Analysis of Private Aircraft Antenna Patterns", 
Report 710964-7, October 1980, The Ohio State University 
ElectroScience Laboratory, Department of Electrical Engineering; 
prepared under Grant NSG 1498 for National Aeronautics and Space 
Administration, Langley Research Center, Hampton, Virginia. 

[40] W.D. Burnside and T. Chu, "Airborne Antenna Patterns Code User's 
Manual", Technical Report 711588-2, March 1980, The Ohio State 
University ElectroScience Laboratory, Department of Electrical 
Engineering; prepared under Contract No. N62269-78-C-0379 for 
Naval Air Development Center. 

[41] R.M. Janes, All the World's Aircraft , 1973-1974, McGraw-Hill Co. 

[42] C.L. Yu and W.D. Burnside, "MLS Airborne Antenna Research", 
Semi-annual Report 2902-22, May 1975, The Ohio State University 
ElectroScience Laboratory, Department of Electrical Engineering; 
prepared under Grant NGL36-003-138 for National Aeronautics and 
Space Administration, Langley Research Center, Hampton, Virginia. 

[43] C.L. Yu and W.D. Burnside, "Research on MLS Airborne Antenna", 
Semi-annual Report 2902-23, April 1976, The Ohio State University 
ElectroScience Laboratory, Department of Electrical Engineering; 
prepared under Grant NGL36-003-138 for National Aeronautics and 
Space Administration, Langley Research Center, Hampton, Virginia. 

[44] C.L. Yu and W.D. Burnside, "Volumetric Pattern Analysis of 
Fuselage-Mounted Airborne Antennas", Technical Report 2902-24, 
April 1976, The Ohio State University ElectroScience Laboratory, 
Department of Electrical Engineering; prepared under Grant 
NGL36-003-138 for National Aeronautics and Space Administration, 
Langley Research Center, Hampton, Virginia. 


327 


[45] W.D. Burnside, R.J. Marhefka and E.L. Felton, "Fuselage-Mounted 
Antenna Code - User's Manual", Technical Report 4583-2 (784583), 
July 1977, The Ohio State University ElectroScience Laboratory, 
Department of Electrical Engineering; prepared under Contract 
No. N62269-76-C-0554 for Naval Air Development Center. 

[46] C.A. Lindberg, "A Shallow-Cavity UHF Crossed-Slot Antenna", IEEE 
Trans. AP-17, September 1969, pp. 558-563. 

[47] G.F. Mason, "Radiation Pattern Measurements for F-16 Antenna 

Systems: Volume I - TACAN Antenna System", Report 16PR383, 

September 15, 1976, General Dynamics, Ft. Worth Division. 

[48] H. Dean Cubley and Haynes Ellis, "Antenna Development for the 
Space Shuttle Orbiter Vehicle", IEEE Trans, on Communications, 

Vol . COM-26, No. 11, November 1978, pp. 1713-1722. 

[49] J.H. Richmond, "Radiation and Scattering by Thin-wire Structures 
in the Complex Frequency Domain", Technical Report 2902-10, July 

1973, The Ohio State University ElectroScience Laboratory, 
Department of Electrical Engineering; prepared under Grant NGL 
36-008-138 for National Aeronautics and Space Administration, 
Washington, D.C. 

[50] J.H. Richmond, "Computer Program for Thin-wire Structures in a 
Homogeneous Conducting Medium", Technical Report 2902-12, January 

1974, The Ohio State University ElectroScience Laboratory, 
Department of Electrical Engineering; prepared under Grant NGL 
36-008-138 for National Aeronautics and Space Administration, 
Washington, D.C. 


328 



